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pJL^ ' Abstract. All 63 discrete X-ray sources brighter than Lx=5xl0 36 erg s 1 in any of four XMM-Newton 

observations of the core of M31 were surveyed for time variability over time-scales of seconds to thousands 
of seconds, and for intensity variations between observations. This population is likely to be dominated 
by low mass X-ray binaries (LMXBs). Analysis of the shapes of power density spectra (PDS) of these 
, sources allows us to determine whether they are accreting at a high or low rate; in the observed frequency 

^ ■ range a broken power law (Type A) PDS indicates a low accretion rate while a simple power law (Type 

B) indicates a high accretion rate. We obtained the 0.3-10 keV luminosities of the sources by modelling 
their spectral energy distributions. The luminosity function for Type A PDS appears to consist of two 
populations, which we tentatively classify as neutron star and black hole LMXBs. We find that 46 sources 
are likely X-ray binaries, 13 with possible black hole primaries. Remarkably, 4 out of 11 LMXB black hole 
candidates appear persistently bright over ~20 years; there are no persistently bright black hole LMXBs 
known in our Galaxy. The lightcurves of 55 X-ray sources had a probability of variability >99.9% in at 
least one observation. Also, 57 of the sources show a luminosity variation <;5cr between observations; 7 of 
• these are classed as transients, since they are absent in at least one observation, and vary in luminosity 

by a factor of <;10. Scaling this sample with the known Galactic LMXBs, we find that the M31 sample 
q . has ~50% of the dippers, Z-sources and transients. 
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■ 1. Introduction telescopes is the largest of any X-ray telescope imaging 

above 2 keV. 

At 760 kpc (van den Bergh 2000) the Andromeda Four XMM-Newton observations were made of the 

Galaxy (M31) is the nearest spiral galaxy to our own. CQre of mi between 20 00 June and 2002 January, as 

The X-ray emission of M31 is dominated by point part of & suryey of thc wholc galaxy Thege observa _ 

sources, mostly X-ray binaries (XBs). As low mass tionS; and other M31 fieldS; represent the best oppor _ 

stars have considerably longer evolution times than tunity fo date for studying time varia bility of X-ray 

high mass stars, low mass X-ray binaries (LMXBs) are sources m an external galaxy 

expected to dominate the central region of the galaxy, „ , , v n i j r ■ 

. . I,,, . i. i 1 he observed X-ray sources m the field of view are 

which is dominated by thc old bulge population: high ,., , , . , r r 

, . . , TT1 , „ r „ . ,., , f likely to be a mixture of supernova remnants, tore- 
mass X-ray binaries (HMXBs are more likely to be , , , , , n j 

v ' ground stars, background galaxies and various accre- 

found m star forming regions such as thc spiral arms. ,. , , . PDin , , , , 

° b ^ tion powered binaries, ihe LPlC-pn detector has a 

X-ray binaries exhibit several types of variability 30'x30' held of view (Struder et al. 2001), while the 

on short time-scales (i.e. milliseconds to hours). The tw0 EPIC-MOS detectors have circular helds of view 

study of such variability in X-ray sources in external w j t h a 30' diameter (Turner et al. 2001), and all three 

galaxies has been limited by the sensitivity of previ- ne i ds f v i ew are ^co-aligned. 

ous observatories, but the three large X-ray telescopes xhe intensity var i at ions on short and long time- 
on board XMM-Newton make such a project possible scales in the x _ ray population of M31 provide vital 
for thc first time; the combined effective area of these clues as to the nature of the sourceS5 and must be ac _ 

counted for in combination with the information ob- 

Send offprint requests to: R. Barnard, e-mail: tained from spectral analysis. Hence a survey has been 
r.barnard@open.ac.uk undertaken, in search of such variability in the com- 
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bined lightcurves of the EPIC MOS and PN detectors 
on board XMM-Newton. Variability of X-ray sources 
over 6 month time-scales is discussed by Osborne et al. 
(2001). 

The core of M31 has been observed many times 
with Chandra (e.g. Kong et al. 2002a; Kaaret 2002; 
Williams et al. 2004a). The high spatial resolution 
of the Chandra images makes them a vital compan- 
ion to the XMM-Newton survey; a single point source 
in an XMM-Newton observation may be resolved by 
Chandra into multiple sources (see e.g. Barnard et al. 
2004). By combining results from both observato- 
ries, we can obtain more information than from either 
alone. 

In this paper, we summarise the known properties 
of our X-ray sample, and discuss their power density 
spectra. In Sect. 2 we review the current understand- 
ing of the X-ray spectra of black hole and neutron star 
LMXBs, while in Sect. 3 we discuss the sources of time 
variability we expected to observe in the M31 LMXBs. 
Most important is Sect. 3.3, which discusses stochas- 
tic variability in disc accretion, and is the foundation 
of our work. In Sect. 4 we show how combining the 
energy and timing analysis allows us to further clas- 
sify the sources. In Sect. 5 we survey the literature on 
variability in Galactic LMXBs with known primaries, 
to support the work outlined in Sect. 3.3. We then 
summarise the observations in Sect. 6 and present the 
results in Sect. 7. Finally, we summarise our findings 
in Sect. 8. 

2. The X-ray spectra of LMXBs 

Black hole LMXBs are thought to have five spec- 
tral states, reviewed comprehensively in McClintock 
& Rcmillard (2004) and Zdziarski & Gierlinski (2004). 
In order of increasing luminosity, these are the quies- 
cent, low/hard, intermediate, high/soft, and very high 
states. The quiescent state is thought to be an exten- 
sion of the low/hard state to low luminosities. The 
spectral energy distribution (SED) of the low/hard 
state is well described by a hard cut off power law, with 
a spectral index, T, ^1.8 and a break at 100-200 keV. 
In contrast, the high/soft state has a soft, thermal 
SED, described by a multicoloured disc blackbody. 
The intermediate and very high states are almost in- 
distinguishable apart from their luminosity, and their 
SEDs may be described by a two component model, 
consisting of a thermal component and a soft power 
law component, with T > 2.4. 

The SEDs of LMXBs with neutron star primaries 
are different from those of black hole LMXBs. The 
continuum spectra of the former are described by 
two spatially distinct spectral components: a thermal 
component, and a non-thermal component. The non- 
thermal component arises from unsaturated, inverse- 
Compton scattering of cool photons on hot electrons 
in an accretion disc corona (ADC); however, the ge- 
ometry of these emission regions remains controver- 



sial. Observations favour a point-like blackbody (i.e. 
the neutron star), and an extended ADC with a 
radius of -10,000-500,000 km (see e.g. Church & 
Balucinska-Church 1995; Barnard et al. 2001; Church 
& Balucinska-Church 2004). At low luminosities, the 
SEDs of neutron star LMXBs may be well described 
by a single Comptonized component (White et al. 
1988), like those of the black hole LMXBs. In contrast 
to black hole LMXBS, however, neutron star LMXBs 
are dominated by the non-thermal component even at 
high luminosities (e.g. White et al. 1988; Barnard et al. 
2003a); the blackbody component contributes ^10% 
at low luminosities and up to 30-50% at higher lumi- 
nosities (Church & Balucinska-Church 2001). Hence, 
a thermally dominated SED at high luminosity is the 
signature of a black hole LMXB (Done & Gierlinski 
2003 2004). 

3. Sources of time variability expected in the 
LMXBs in M31 

Long term variability in the point X-ray sources of 
M31 has been well studied in the last 5 years. Two 
ROSAT PSPC catalogues of M31 point sources were 
compiled (Supper et al. 1997 2001, hereafter S97 and 
S01 respectively), with a combined catalogue of 560 
point sources. The 1991 survey revealed 15 variable 
sources and 9 transient sources when compared with 
a previous Einstein catalogue (Trinchieri & Fabbiano 
1991), while comparison between the 1992 and 1991 
surveys revealed 34 long term variables and 8 pos- 
sible transients (Supper et al. 2001). The 1999-2001 
Chandra observations reveal variability from obser- 
vation to observation in —50% of the X-ray sources 
(Kong et al. 2002a, hereafter known as K02), while 
Trudolyubov & Priedhorsky (2004) report variabil- 
ity in 80% of globular cluster sources over 10 XMM- 
Newton and ~40 Chandra observations. 

Galactic LMXBs can exhibit a variety of phe- 
nomena resulting in luminosity and colour varia- 
tions over timescales of milliseconds to hours. Two 
of these phenomena, namely dipping and Z-source 
movement, should be identifiable in the lightcurves of 
M31 LMXBs. These will be discussed in Sect. 3.1 and 
Sect. 3.2 respectively. Additionally, the disc accretion 
process results in stochastic variability, discussed in 
Sect. 3.3. 

3.1. Periodic intensity dips due to photo-electric 
absorption 

Several high inclination Galactic LMXBs experience 
periodic modulation of their lightcurves and SEDs, 
due to obscuration of the X-ray source by mate- 
rial that is raised out of the body of the disc (see 
e.g. White & Swank 1982; Church & Balucinska- 
Church 1995). Photoelectric absorption leads to pref- 
erential removal of soft photons, hence hardening the 
SED. Dipping has been reported in XMM-Newton 
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observations of two LMXBs in the core of M31: 
XMMUJ004314.1+410724 (Bo 158, Trudolyubov 
ct al. 2002) and XMMU J004308. 6+411247 (Mangano 
et al. 2004), on periods of 2.78 and 1.78 hr respec- 
tively. The dipping behaviour of Bo 158 is particu- 
larly interesting, as the depth of dipping is ~82% in 
the 2002, January, observation, yet the dip modulation 
is SS10% in the 2001, June observation (Trudolyubov 
et al. 2002). In Barnard et al. (2006), we show that 
the evolution in the dipping behaviour may be due 
to warped disc precession; we modelled the disc with 
three dimensional smoothed-particle hydrodynamics, 
and find that the disc is indeed expected to precess. 

3.2. Z-source movement 

The Galactic neutron star LMXBs are divided into Z- 
sources and atoll sources, based on their colour-colour 
diagrams, or colour-intensity diagrams, and timing 
properties (Hasinger & van der Klis 1989). The six 
Galactic Z-sources all have luminosities in excess of 
10 38 erg s _1 , and trace out a three-branched "Z" in 
their colour-colour diagrams. The three branches are 
called the horizontal branch, normal branch and flar- 
ing branch, and all Galactic Z-sources exhibit a ~6 
Hz quasi-periodic oscillation (QPO) on the normal 
branch (Hasinger & van der Klis 1989). Branch move- 
ment of Z-sources results in significant changes in the 
SED and luminosity over time-scales of minutes to 
hours (e.g. Barnard et al. 2003a). Atoll sources have 
now been shown to also exhibit Z-shaped colour- 
colour diagrams (Gierlihski & Done 2002). However, 
atoll sources trace their Z-tracks over timescales of 
~months, requiring intensity variations by a factor of 
~80, whereas Z-sources only vary in intensity by a 
factor of ~2-3, and traverse the Z over a few days 
(Muno ct al. 2002). The difference between Z-sources 
and atoll sources is therefore real, but of unknown 
origin. While most LMXBs have orbital periods ^10 
hr, the Z-sources have orbital periods ^20 hr (e.g. 
Liu et al. 2001); hence, the secondaries in Z-source 
are likely to be larger and/or more evolved. It is be- 
lieved that the difference between Z-sources and atoll 
sources may be due to some vital difference in their 
secondaries. 

In Barnard et al. (2003c), we identified a likely 
Z-source in M31, RX J0042. 6+4115. We produced a 
colour-intensity plot for four XMM-Newton observa- 
tions of the source, defining the colour as the inten- 
sity ratio of the 4.0-10 kcV and 0.3-2.0 keV bands, 
and found that the data clustered in three regions, 
roughly describing a Z. The 0.3-10 keV SED of 
RX J0042. 6+4115 required a two component spectral 
model, consisting of a blackbody component and a 
power law component arising from Compton scat- 
tering, with an unabsorbed luminosity of ^5xl0 38 
erg s" 1 in the 0.3-10 keV band; RX J0042. 6+4115 
is a good M31 analogue for the Galactic Z-sources 
(Barnard et al. 2003c). 



3.3. Stochastic variability in disc accretion 

Van der Klis (1994) showed that the power density 
spectra (PDS) of low mass X-ray binaries (LMXBs) 
depend more on the accretion rate than on the na- 
ture of the primary. At low accretion rates, the frac- 
tional rms variability is high (a few tens of percent, 
van der Klis 1995) and the PDS are well described 
by a broken power law that changes in spectral in- 
dex, 7, from ^0 to ~1 at frequencies higher than a 
certain break frequency; the break occurs at 0.01-1 
Hz (van der Klis 1994). At higher accretion rates, the 
rms variability is only a few percent, and the PDS are 
characterised by a power law with 7 ~1-1.5. Van der 
Klis (1994) proposed that LMXBs switch from low to 
high accretion rate behaviour at some constant frac- 
tion of the Eddington limit, and suggested a transition 
at 1% Eddington. Since we cannot observe the accre- 
tion rate, M, directly, we must use the luminosity to 
trace the evolution of the PDS with M. While the lu- 
minosity is a good indicator of M, the correspondence 
is not simple (see e.g. Done & Gierlihski 2004). We 
dehne l c as 



where L c is the luminosity of transition and L^dd is 
the Eddington limit. 

We denote the low accretion rate variability (bro- 
ken power law PDS) as Type A and high accretion 
rate variability (simple power law PDS) as Type B 
(Barnard, Kolb, & Osborne 2004). We realised that 
if the transition from Type A to Type B occurred at 
some constant l c , then black hole LMXBs would ex- 
hibit Type A variability at higher luminosities than 
neutron star LMXBs (Barnard et al. 2003b 2004). We 
have some evidence that l c ~0.1 in the 0.3-10 keV 
band, from Galactic LMXBs with known primaries, 
and the globular cluster X-ray sources in M31. We 
present results from the literature in Sect. 5, and evi- 
dence from the M31 globular cluster X-ray population 
in Sect. 7.2. 

4. Combining SED and Timing analysis 

McClintock & Remillard (2004) recently reviewed the 
spectral and timing properties of Galactic black hole 
binaries in their different states. Ten years after the 
van der Klis (1994) paper, the general scheme of black 
hole behaviour still holds. For the hard (low/hard) 
state, McClintock & Remillard (2004) observe the 
SED to be dominated by a power law component with 
T -1.5-2.1, while the PDS is of Type A, with r.m.s. 
variability ~0. 1-0.3 in the 0.1-10 Hz band. For the 
thermally dominated (high/soft) state, a disc black 
body with temperature 0.7-2 keV contributes ^90% 
of the flux in the 0.3-10 keV band; Type B PDS are ob- 
served with fractional r.m.s. variabilities of ^0.01-0.06 
over the 0.1-10 Hz frequency range. McClintock & 
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Table 1. Classifications of the X-ray population of M31 by power density spectrum (PDS) and spectral energy distri- 
bution (SED). The SED is classified by the photon index of the power law component, T, or by the temperature of the 
thermal component, YTbb- 



PDS 


SED 


Classification 


Type A 


r 5-2.1 


NS or BH LMXB in low/hard state 


Type A 


r ~2.4-3.0 


BH LMXB in very high/ 






intermediate (steep power law) state 


TypeB 


BB dominated 


BH LMXB in high/soft state 




kT BB ~0.7-2 keV 




TypeB 


Not BB dominated 


NS LMXB in high state 



Rcmillard (2004) replace "very high state" and "inter- 
mediate state" with the "steep power law" state. The 
SED contains a power law component with T ~2. 4-3.0 
which contributes >20% of the 2-20 keV emission. If 
the power law component contributes 20-50%, then a 
QPO is seen in the PDS at ^0.1-30 Hz; however, if 
the power law component contributes >50%, then no 
QPO is seen. The PDS are of Type A with break fre- 
quencies ~0.1-10 Hz, i.e. at higher frequencies than 
in the low hard state. There appears to be no con- 
sensus on the r.m.s. variability of the PDS associated 
with the very high/intermediate state; however, van 
der Klis (1995) states that the r.m.s. variability lies 
between those of the low/hard and high/soft states. 

Joint analysis of the SEDs and PDS of the M31 
X-ray sources therefore allows greatly improved clas- 
sification of the population. In Table 1, we summarise 
our classification scheme. 

5. PDS transitions in Galactic LMXBs with 
known primaries 

We have performed a literature survey of Galactic 
LMXBs, with the goal of constraining l c . We consid- 
ered LMXBs with known primaries that have either 
shown transitions from Type A to Type B, exhib- 
ited Types A and B within a narrow luminosity range, 
or have peaked during outburst without making the 
transition to Type B. The motivation is to discover 
whether the transitions occur at some constant frac- 
tion of the Eddington limit; if true, then the transition 
luminosity is a probe of the primary mass. Since we 
are unable to measure the bolometric luminosities of 
M31 X-ray sources, we have converted all luminosi- 
ties presented in the literature to luminosities in the 
0.3-10 keV band. 

5.1. l c for a Galactic neutron star LMXB 

4U 1705-44 is a Galactic LMXB that exhibits X-ray 
bursts, and hence contains a neutron star (Langmeier 
et al. 1987). It exhibited a Type A PDS in the faintest 
of four EXOSAT observations, and a Type B PDS in 
the next faintest; the respective 1-20 keV fluxes were 
~T.3xl0~ 9 and ~T.8xlO~ 9 erg cm~ 2 s _1 (Langmeier 
et al. 1987 1989). We obtained 0.3-10 keV fluxes for 



these observations using the best fit spectral models 
obtained by Langmeier et al. (1987), yielding fluxes of 
1.2 and 1.9 xl0~ 9 erg cm~ 2 s^ 1 respectively. Hence, 
an accurate distance could yield a tight constraint 
on l c . The distance to 4U 1705—44 has been estimated 
using X-ray bursts as standard candles (see Kuulkcrs 
et al. 2003); Christian & Swank (1997) find a distance 
of 11 kpc from Einstein data, while Cornelisse et al. 
(2002) obtain a distance of 8.6 kpc using BeppoSAX. 
If we assume that the distance lies between these two 
values, then l c is consistent with 0.10±0.04. 

5.2. l c for Galactic black hole LMXBs 

Of the 18 confirmed Galactic black hole X-ray bina- 
ries, 3 are high mass X-ray binaries (HMXBs), and 
are persistently bright, and the rest are transient 
LMXBs. In addition, there are a further 22 transient 
LMXBs that are black hole candidates (McClintock & 
Remillard 2004). These 37 LMXB black hole or black 
hole candidate transients make up ^50% of the known 
Galactic LMXB transients (see e.g. in't Zand et al. 
2004). In general outbursts last several months and 
the X-ray luminosity can increase by a factor of 10 7 ; 
outbursts are on average separated by years of quies- 
cence (Chen et al. 1997; in't Zand et al. 2004). 

Transient LMXBs are hysteretic, in that the hard 
to soft transition during the rise of the outburst oc- 
curs at a higher luminosity than the soft to hard tran- 
sition during decay (Miyamoto et al. 1995; Maccarone 
& Coppi 2003). We were interested in obtaining the 
maximum value for / c , as this would be most sensitive 
to the Eddington limit. Hence, estimates of l c in black 
hole LMXBs were restricted to those that had either 
been observed during the rise of the outburst, or ex- 
hibited outbursts where the transition from low/hard 
state to high/soft state was not made (see Brocksopp 
et al. 2004, for a review). Unfortunately, most X- 
ray observations of Galactic black hole LMXBs during 
outburst have been during the decay phase, where the 
spectral transition occurs at I < l c due to hysteresis, 
and are hence unsuitable for this work. 

Nine black hole LMXBs have exhibited outbursts 
where they remained in the low/hard state (for a re- 
view see Brocksopp et al. 2004). Five of these have 
published Type A PDS, distances and mass estimates; 
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Name 


L A a 


Mass / M 


Min l c h 


GS 2023+338 


>30 (2-37 keV) c 


~12 d 


0.06±0.01 e 


GRSJ1719-24 


3 (3-200 keV) / 


~5 / 


0.06+0.01 9 


GRS 1737-31 


3 (2-200 keV) h 


4-14 


0.015+0.005* 


GRO 0422+32 


6.4+0.8 (3-200 keV) j 


4.0±1.0 fc 


0.10+0.03' 


XTE J1118+480 


0.04 (2-12 keV) 1 " 


6.0-7.7" 


0.00067+0.00018° 



a Highest luminosity where Type A PDS is observed / 10 37 erg s -1 ; "minimum l c in 0.3-10 keV band; c (Miyamoto et al. 
1992; Oosterbroek et al. 1997, and references within); d Shahbaz et al. (1994); e (assuming power law X-ray spectrum 
with spectral index a = 1—1.4 Miyamoto et al. 1992); ^Revnivtsev et al. (1998), assuming distance of 2.8 kpc (della 
Valle et al. 1994); 9 a = 2-2.3 (van der Hooft et al. 1996); ''Cui et al. (1997); l a ~1.7 (Cui et al. 1997); 'obtained 
by comparison of 40-150 keV flux with GRS 1737—31, which has a similar a (van der Hooft et al. 1999) and distance 
(Gelino & Harrison 2003); fc Gelino & Harrison (2003); l a -2.1 (van der Hooft et al. 1999); m Flux = 42 mCrab (Hynes 
et al. 2003), distance^!. 71+0.05 kpc (Chaty et al. 2003); "Wagner et al. (2001); °a = 1.78 (McClintock et al. 2001) 



we have used published results to obtain the corre- 
sponding minimum value of l c in the 0.3-10 keV band. 
The results are presented in Table 2; the mass of 
GRS 1737—31 is given by the mass range of known 
black holes. The minimum l c values in Table 2 cover 
the range 7xl0~ 4 -0.10, with l c >0.06 for 3 out of 5 
sources. It is reasonable to conclude that l c ^0.05 in 
the 0.3-10 keV band. 

Two of the confirmed black hole LMXBs, 
GX 339-4 (Zdziarski et al. 2004) and 
XTE J1550-564 (Rodriguez et al. 2003), have recently 
been caught during the rise of the outburst by the 
RXTE-ASM, allowing monitoring of the entire out- 
burst by the main instruments of RXTE. Zdziarski 
et al. (2004) present results from two outbursts of 
GX 339-4 that were observed with RXTE, and a 
lightcurve spanning 1987-2004, covering —15 out- 
bursts. They found that the luminosity of the tran- 
sition between spectral states during the rise of the 
outburst depended on the history of the disc: the tran- 
sition during the first outburst observed by RXTE was 
at a bolometric luminosity, L bol , of 7% Eddington, 
while the transition during the rise of the second out- 
burst was at Lboi ~ 20% Eddington; these correspond 
to 10-3-wkeV _ 0.04 and 0.10 respectively. Prior to 
the first outburst, GX339— 4 was quasi-persistently 
bright with frequent intensity dips, whereas the sec- 
ond outburst was preceded by a —1000 day low state 
(Zdziarski et al. 2004). The second outburst is more 
representative of a canonical outburst, and has simi- 
lar characteristics to the April- June 2000 outburst of 
XTE J1550— 564, where the spectral state transition 
during the rise also occurred at Lboi ~20% Eddington, 
corresponding to —10% in the 0.3-10 keV band. The 
spectral transition was not observed in its 1998 out- 
burst, instead the very high, high and intermediate 
states were observed (Sobczak et al. 1999). We note, 
however, that Zdziarski et al. (2004) and Rodriguez 
et al. (2003) identified the transitions from the chang- 
ing SEDs, not their PDS. We assume that the tran- 
sition from Type A to Type B variability is simulta- 



neous with the transition from the low/ hard to the 
high/soft state, although this may not be the case. 

5.3. An empirical value for l c 

Our considerations show that eight Galactic LMXBs 
with known primaries (7 black holes, 1 neutron star) 
are consistent with l c —0.1 in the 0.3-10 keV band. It 
is beyond the scope of this paper to conduct a compre- 
hensive survey of PDS in Galactic LMXBs. Instead, 
we intend these results to be a useful starting point, 
prior to a thorough investigation of PDS transitions 
in Galactic X-ray binaries with known primaries. 

Esin et al. (1997) modelled the hard to soft tran- 
sition using observations of the black hole transient 
Nova Muscae 1991. They found that the transition 
occurs at a critical accretion rate as a fraction of the 
Eddington limit, m cr i t that depends on the disc vis- 
cosity a: m C rit ~1.3 a 2 . They assume a = 0.25, giving 
TOcrit ^0.08. Esin et al. (1998) then successfully ap- 
plied the model to CygnusX-1, GROJ0422+32 and 
GROJ1719— 24, obtaining m cr i t —0.1. However, it is 
uncertain whether l c = m cr it- 

We note that Esin et al. (1998) find that the SED 
of the GRO 0422+32 at the peak of its outburst is 
consistent with an accretion rate of — m cr j t ; this is in- 
teresting because we find the 0.3-10 keV luminosity of 
GRO 0422+32 at that time to be 10+3% Eddington. 
Results from the other six black hole LMXBs are also 
consistent with l c —0.1 in the 0.3-10 keV band, al- 
though they cannot constrain l c . For the remainder of 
the paper, we assume that l c = 0.1 for all LMXBs. 

The maximum mass for a neutron star is generally 
accepted to be 3.1 M (e.g. Krishan & Kumar 1978), 
but could be as high as 6 M 0) depending on the equa- 
tion of state (Srinivasan 2002). We must appeal to ob- 
servations: all measured neutron star masses are less 
than 2 M , and most are consistent with 1.35+0.04 
Mq (Srinivasan 2002, and references within). We will 
assume a maximum neutron star mass of 3.1 Mq. This 
allows us to postulate a maximum luminosity for Type 
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Table 3. Journal of XMM-Newton observations of the 
M31 core. The observation number, date, total exposure 
(Exp), exposure of the good time interval (GT) and filter 
are shown. 



Observation 


Date 


Exp 


GT 


Filter 


1 


25/07/00 


34 ks 


27 ks 


Medium 


2 


27/12/00 


13 ks 


13 ks 


Medium 


3 


29/06/01 


56 ks 


25 ks 


Medium 


4 


06/01/02 


64 ks 


64 ks 


Thin 



A variability in a neutron star LMXB of ^4x10 erg 
s- 1 in the 0.3 10 kcV band. 

6. The Observations 

Four observations of the core of M31 were made with 
XMM-Newton, separated by 6 months; a journal of 
the observations is presented in Table 3. Observations 
1 and 3 suffered intervals of background flaring; for 
each of these observations, a single flare- free good time 
interval was used. While previous surveys used many 
short pointed observations, the four XMM-Newton ob- 
servations were long and uninterrupted; together with 
the unprecedented effective area of XMM, the obser- 
vations yielded up to 40 times the counts of the best 
previous observations. 

Data reduction was performed using the XMM- 
Newton SAS version 6.0.0, and the products were 
analysed with XANADU and the FTOOLS software 
suite. We conducted a timing survey of the 63 bright- 
est X-ray sources in these four XMM-Newton observa- 
tions of the central region of M31, extracting the PDS 
and 0.3-10 keV luminosity of each source, assuming a 
distance to M31 of 760 kpc (van den Bergh 2000). 

We selected our sources by intensity; the mean in- 
tensity of the background subtracted, 0.3-10 keV pn 
lightcurve had to exceed 0.02 count s _1 in at least one 
observation; this ensured at least ~250 source counts 
in the pn spectrum from the shortest observation. For 
each source, data were extracted from a circular region 
that was centred on the source; those sources within 
8' of the nucleus had extraction regions with 20" radii, 
due to crowding. Sources further than 8' from the nu- 
cleus had extraction regions of 40". A background re- 
gion was chosen for each source to be a source-free 
region with the same size as the source region, on the 
same CCD and at a similar off-axis angle. 

Source and background lightcurves with 2.6 s time 
resolution were obtained for each source from the 
three EPIC detectors: MOS1, MOS2 and pn. The 
background-subtracted lightcurves were corrected for 
vignetting and summed, giving a combined EPIC 
lightcurve for every source, for each of the four 
XMM-Newton observations. PDS were then obtained 
for each lightcurve using the powspec program in 
FTOOLS, averaging over intervals of 333 s and 666 
s with time bins of 5.2 s (intervals of 64 and 128 bins 



respectively). PDS were also made for the combined 
EPIC background lightcurves. This allowed us to de- 
termine whether the variability observed in the source 
PDS was caused by variation in the background. 

The output PDS were then fitted with power law 
and broken power law models, using custom software. 
The power law model was parameterised by the spec- 
tral index, 7, and normalisation, while the broken 
power law model was characterised by the spectral 
index changing from a to /? at the break frequency 
v c , normalised by the power at v c . To ensure that the 
broken power law fits described Type A PDS, we set 
a = 0. We then conducted an F-test on the best fit 
power law and broken power law models, to determine 
whether the broken power law gave a significantly bet- 
ter fit. 

Those PDS that were well fitted by broken power 
law models, when power law models could be rejected 
at a level J;99.9%, were classified as Type A. We also 
classified as Type A those PDS where the probabil- 
ity of improved fitting with a broken power law was 
>99.9%. PDS that were well described by a simple 
power law, and clearly displayed power in excess of 
the Poisson noise, were classified as Type B. Those 
PDS which exhibited no excess power were classed as 
flat. Any PDS that were not classifiable by any of the 
above means were rejected. We only accepted PDS 
classifications where the 64 bin and 128 bin PDS of 
the same lightcurve were classified as the same type 
(i.e. both Type A, both Type B or both flat). 

We simulated random lightcurves with Type B 
PDS following the method of Timmcr & Kocnig (1995) 
and normalised the intensities to 0.5 count s _1 (a typ- 
ical combined EPIC intensity of a bright X-ray binary 
in M31). We found that the power often dipped be- 
low the Poisson noise level in the observed frequency 
range. Hence, a flat PDS with an average source inten- 
sity that was higher than the faintest Type A source 
in the same observation was consistent with Type B 
but not Type A. 

Source and background SEDs were extracted for 
each EPIC-pn observation of every source in 4096 
channels, with 5 eV binning. A response matrix was 
generated for each source, and corresponding ancil- 
lary response files were also made. Counts outside the 
0.3-10 keV range were rejected, and the spectra were 
grouped according to the source intensity: spectra ex- 
ceeding 500 counts were grouped for a minimum of 50 
counts per bin, and spectra containing less than 500 
counts were grouped to a minimum of 20 counts per 
bin. 

7. Results 

The 63 sources that satisfied the intensity selection in 
at least one of the four observations are catalogued 
in Table 4. The J2000 coordinates are given, along 
with their designations in the first and second ROSAT 
PSPC catalogues (S97, S01) and the Chandra cata- 
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Table 4. Coordinates, corresponding catalogue number in S97, SOI, K02 and W04 catalogues, along with properties 
recorded in these surveys, PDS types observed and the most likely primary when known. The codes in the properties 
column are explained in the text. References for previously known properties are given in the foot notes; a property 
without a reference is first reported in this work. 



Qni 1 roc* 


R.A. 


Dec. 


S97 


S01 


K02 /W04 


vnmavt" i pq 


PDS 


_L 1 1111CL1 y 


si 


00:42:05.71 


+41:13:29.7- 7 ' 






r3-125 


T c v 


A - 




S2 


00:42:07.619 


+41:18:15. 17 d 


141 


139 


r3-61 


t a 


- 




S3 


00:42:08.952 


+41:20:48.42 d 


142 


142 


r3-60 


V 


A - 




S4 


00:42:09.372 


+41:17:45.63 d 






r3-59 


glob d t l v 


B - 




S5 


00:42:12.026 


+41:17:58.86 d 






r3-54 


t c glob d v 






S6 


00:42:13.017 


+41:18:36.73 d 


146 


144 


r3-52 


t a fg a v 


- 




S7 


00:42:15.011 


+41:12:34. 23 d 






r3-50 


t d V 


A B 


bh 


S8 


00:42:15.571 


+41:17:21. ll d 


151 




r3-47 


e v 


A 




S9 


00:42:15.86 


+41:01:14.7' 


150 


147 


sl-7 


t a glob" 


- 




S10 


00:42:18.241 


+41:12:23. 53 d 


155 


153 


r3-45 


t a V 


- 




Sll 


00:42:18.534 


+41:14:01.69 d 


158 


154 


r3-44 


t a glob a v 


- 




S12 


00:42:21.382 


+41:16:01.32 d 






r3-42 


t d v 


A 


bh 


S13 


00:42:22.316 


+41:13:33.99 d 






r3-40 


t d V 


A 




S14 


00:42:22.841 


+41:15:35. 14 d 


163 


159 


r3-39 


t d v e 


A 9 


bh 9 


S15 


00:42:25.941 


+41:19:15. 27 d 






r2-36 


glob d t d v 


B 




S16 


00:42:26.047 


+41:25:52. 74 d 


167 


166 


r3-87 


V 


A B 




S17 


00:42:28.089 


+41:09:59.84 d 






r3-36 


t d V 


A - 




S18 


00:42:28.190 


+41:12:22. 76 d 


172 


169 


r2-35 


e v 






S19 


00:42:28.789 


+41:04:34.98 d 


173 


170 


r3-lll 


e v 


B 




S20 


00:42:31.045 


+41:16:21. 74 d 


176 


174 


r2-34 


t d 


- 




S21 


00:42:31.154 


+41:19:39. 19 d 


175 


175 


r2-33 


glob a t d v 


A - 




S22 


00:42:31.979 


+41:13:14. 24 d 


177 


177 


r2-32 


t d m(2) v 


B 




b23 


1)0:42:34.361 


i a i . 1 o .nn add 
+41:18:09.60 






~o on 
rz-z9 


1 V 


A 




S24 


00:42:35.121 


+41:20:06.09 d 


182 




r2-27 


t d V 


A - 




S25 


00:42:38.503 


+41:16:03.80 d 


184 


184 


r2-26 


t d zs h V 


B - 


ns 


S26 


00:42:39.451 


+41:14:28. 52 d 






r2-25 


t d v 






S27 


00:42:40.121 


+41:18:45.38 d 






r2-24 


t d v 


A 




S28 


00:42:41.566 


+41:21:06.02 d 






r3-31 


t d v 


A 




S29 


00:42:43.225 


+41:13:19.48 d 






r2-19 


V 


A - 




S30 


00:42:44.766 


+41:ll:37.76 d 


195 


194 


r3-29 


t d v 


B 


ns 


S31 


00:42:47.089 


+41:16:28. 65 d 


198 


195 


rl-2 


t c pn d m(4) v 


A 




S32 


00:42:48.450 


+41:25:23. 10 d 


201 


198 


r3-25 


t d v 


B - 





a S97, S01; b Kahabka (1999); c Osborne et al. (2001); d K02; e Kaaret (2002); / Trudolyubov et al. (2002); 9 Barnard 
ct al. (2003b); ^Barnard, Kolb, & Osborne (2003c); ^Mangano et al. (2004); J W04; fc Barnard, Kolb, & Osborne (2004); 
'Trudolyubov & Priedhorsky (2004) 



logucs of K02 and Williams et al. (2004a), hereafter 
referred to as W04; known properties of each source 
are given, with separate columns for the observed PDS 
types and the most likely primaries where known. The 
coordinates are taken from K02, or W04 for sources 
not in the K02 catalogue. Intensity variability over 
short ('-TOO s) timescales is indicated by v, and X-ray 
sources that vary in luminosity by more than 5a be- 
tween observations are labelled t. We define transients 
(labelled T) as sources that vary in intensity by a fac- 
tor of ^10, and are absent in at least one observation. 
Foreground sources are indicated by fg, globular clus- 
ter sources by glob, supernova remnants by snr and 
supersoft sources by sss. Dipping is indicated by dip, 
Z-source branch movement by zs, pulsations by pulse, 
and positional coincidence with planetary nebulae by 
pn. A likely neutron star primary, as prescribed in 



Table 1, is indicated by ns, while a black hole can- 
didate is indicated by bh. Radio counterparts arc in- 
dicated by r, and are taken from the catalogues of 
Walterbos et al. (1985) and Walterbos et al. (2002). 
K02 found positional coincidences with the globular 
clusters using the catalogues of Battistini et al. (1987), 
Magnier (1993), and Barmby (2000). The planetary 
nebulae were assigned using the catalogues of Ford 
& Jacoby (1978), and Ciardullo et al. (1989) — K02 
comment that these X-ray sources are 3 orders of mag- 
nitude brighter than for PN. Williams et al. (2004b) 
obtained the X-ray positions of the six known PN asso- 
ciations from a 37 ks Chandra ACIS observation; they 
found that the planetary nebulae are not responsible 
for the X-rays, but are likely to be associated with 
the X-ray emitter. The reasons for these associations 
are not clear. The SNR were identified from d'Odorico 
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Table 4. continued 



Source 


R.A. 


Dec. 


S97 


S01 


K02/W04 


Properties 


PDS 


Primary 


S33 


00:42:48.52 


+41:15:21.4 J ' 


200 


197 


rl-1 


t m(2) v 


A B — 




S34 


00:42:52.450 


+41:18:54. 75 d 


206 


200 


r2-13 


t d v e 


B - 




S35 


00:42:52.450 


+41:15:40. 20 d 


208 


203 


r2-12 


t a sss b V 


- 




S36 


00:42:54.3 


+41:30:49 










- 




S37 


00:42:54.859 


+41:16:03.46 d 


210 




r2-ll 


v t d r? 


A - 


bh 


S38 


00:42:55.316 


+41:25:56.60 d 


211 


206 


r3-23 


snr a t v 


A 




S39 


00:42:55.550 


+41:18:35.44 d 






r2-9 


t d m(2) 


A - 




S40 


00:42:57.854 


+41:11:04. 59 d 


213 


209 


r3-22 


V 


A - 


bh 


S41 


00:42:58.257 


+41:15:29.46 d 






r2-7 


t d V 


A 




S42 


00:42:59.594 


+41:19:19. 72 d 


217 


210 


r2-6 


glob a t d - 






S43 


00:42:59.803 


+41:16:06.01 d 






r2-5 


t d glob d v 


A - 




S44 


00:43:00.87 


+41:30:07.7 J 


218 


212 


nl-76 


glob a 


- 




S45 


00:43:02.958 


+41:20:42. 54 d 






r3-21 


pn d t v 


A B 




S46 


00:43:03.089 


+41:10:15. 18 d 


211 




r3-20 


t V 


A B 




S47 


00:43:03.163 


+41:15:28.00 d 


226 


214 


r2-3 


glob a T d m(2) fc 


A k _ 


bh k 


S48 


00:43:03.231 


+41:21:22.42 d 


222 


216 


r3-19 


glob a t d v 


A 




S49 


00:43:03.812 


+41:18:05. 23 d 


223 


217 


r2-2 


glob a t d v 


A 




S50 


00:43:04.186 


+41:16:01.62 d 






r2-l 


t d V 


A B 




S51 


00:43:05.66 


+41:17:03.3 J 






r2-69 


T j v 


A 


bh 3 


S52 


00:43:08.63 


+41:12:50. V 






r3-17 


v l dip 1 t 


A B - 




S53 


00:43:09.791 


+41:19:01. 22 d 


226 


222 


r3-16 


t a T d rv 


A 


bh 


S54 


00:43:10.587 


+41:14:51.55 d 


228 


223 


r3-15 


glob a t d v 


- 




S55 


00:43:14.245 


+41:07:25.42 d 


229 


227 


r3-112 


glob" dip 1 t f v 


B 




S56 


00:43:18.773 


+41:20:18. 52 d 


235 


235 


r3-8 


t a sss c v e 


A 




DO i 


uu.4o.iy.oz 


\A 1.17. Cft 7J 

-T'kL.l ( .DO. ( 








sss puise l v Tnyz) 


A 




S58 


00:43:27.763 


+41:18:29.83 d 


240 


241 


r3-63 


snr a t a sss e v 


A 




S59 


00:43:29.038 


+41:07:49. ll d 


241 


242 


r3-103 


t V 


A - 


bh 


S60 


00:43:32.382 


+41:10:41.66 d 


243 


244 


r3-3 


t V 


A 


bh 


S61 


00:43:34.332 


+41:13:23. 72 d 


244 


247 


r3-2 


t V 


A 




S62 


00:43:37.191 


+41:14:407 d 


247 


249 


r3-l 


glob a t v 






S63 


00:43:54.2 


+41:16:44 c 


256 


262 


nl-79 


t V 


A 


bh 



a S97, S01; 6 Kahabka (1999); c Osborne et al. (2001); d K02; e Kaaret (2002); / Trudolyubov et al. (2002); 9 Barnard et al. 
(2003b); ^Barnard et al. (2003c); l Mangano et al. (2004); j W04; fc Barnard et al. (2004); 'Trudolyubov & Priedhorsky 
(2004) 



ct al. (1980), Braun & Walterbos (1993), and Magnier 
et al. (1995). Point sources in XMM-Newton that are 
resolved by Chandra into multiple sources are desig- 
nated m{x) 1 where x is the number of sources observed 
by Chandra. Finally, the types of variability are indi- 
cated: A for Type A, B for Type B, and — for fiat 
PDS. If more than one type of variability is observed 
from a given source, then they are observed in different 
observations. 



In the following sections, we first describe the prop- 
erties of PDS observed from the population as a whole 
(Sect. 7.1). We present luminosity functions of Type 
A, Type B and flat PDS, and use them to iden- 
tify possible black hole candidates. We then discuss 
new insights into the globular cluster X-ray sources in 
Sect. 7.2, supernova remnants in Sect. 7.3, our black 
hole candidates in Sect. 7.4, and supersoft sources in 
Sect. 7.5. Our results shed no new light on the known 
dipping sources. 



7.1. Power density spectra 

We provide a description of the PDS and SEDs for 
each observation of every source in Appendix A. We 
were able to classify 107 of the 252 power density spec- 
tra. Examples of Type A and Type B PDS are pro- 
vided in Fig. 1 and Fig. 2 respectively. The axes of 
these figures are logarithmic, and the power is Leahy 
normalised so that the Poisson noise has a power value 
of 2. The power is measured in units of count s _1 
Hz- 1 . 

In Fig. 1 we present a 128 bin PDS of S45 from 
Observation 4, averaged over 96 intervals. Best fit 
power law (dashed line) and broken power law (solid 
line) models are presented, with x 2 /dof = 109/12 and 
3/11 respectively. Figure 2 shows the equivalent PDS 
for S32 in Observation 4. The PDS is well described by 
a simple power law, with x 2 /dof = 11/12. The spec- 
tral index, 7 = —0.07+0.01, is shallower than expected 
for Galactic LMXBs in a Type B state; however, it 
is consistent with the PDS of our simulated Type B 
lightcurves. 
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TYPE A 




2x10 



5xl0~° 0.01 0.02 
Frequency (Hz) 



0.05 



Fig. 1. Example Type A PDS, exhibited by S45 in 
Observation 4; the PDS is averaged over 96 intervals of 
128 bins, and grouped. The axes are log-scaled. The PDS 
is Leahy normalised, so that the Poisson noise has a power 
of 2. Fits of power law (PL) and broken power law (BPL) 
models to the PDS are shown; the power law fit is clearly 
unacceptable, while the broken power law fits well. 

TYPE B 



^ in 

a 



i 

o 
Q. 



PL: xVdof = 11/12 




2x10 



5xlO _-J 0.01 0.02 
Frequency (Hz) 



0.05 



Fig. 2. Example Type B PDS exhibited by S32 in 
Observation 4; the PDS is averaged over 96 intervals of 
128 bins, and grouped. The axes are log-scaled. The PDS 
is Leahy normalised, so that the Poisson noise has a power 
of 2. 



7.1.1. Type A PDS 

We observed 53 Type A PDS in total, over a lumi- 
nosity range of ~5-400x 10 36 erg s _1 in the 0.3-10 
keV band. The break frequencies of our Type A PDS 
range over ~0. 02-0. 06 Hz, which is in the range 0.01-1 
Hz found in Galactic LMXBs by van der Klis (1994). 
However, we note that these break frequencies are 
systematically lower than those observed in Galactic 
LMXBs, which are typically ~0.1 Hz and usually in 
the range 0.04-0.4 Hz (van der Klis 1995). This may 
suggest that the M31 X-ray population may differ sys- 
tematically from the Galactic LMXB population, al- 
though we do not believe this to be the case neces- 
sarily. The observed range of (3 was ~0.2-1.9; however 



5S 




36.5 



37 37.5 38 

Log (Luminosity / erg s _1 ) 



38.5 



Fig. 3. Cumulative luminosity functions for Type A (grey) 
and Type B/flat PDS (black) in the 0.3-10 keV band. 
Dotted lines indicate luminosity functions where each 
source is only counted once, while solid lines indicate 
that all the data are used. The dotted vertical line indi- 
cates a luminosity of 4xl0 37 erg s _1 , the theoretical max- 
imum luminosity for Type A variability in a neutron star 
LMXB(see text). Indeed, we see that the LF above 4x 10 37 
erg s~ is considerably steeper than below. The shallow 
and steep parts of the LF could represent the neutron star 
and black hole populations, respectively. 



ao 

o 
j 




36.5 



37 37.5 38 

Log (Luminosity / erg s _1 ) 



Fig. 4. Cumulative luminosity functions for Type 
in black hole candidates (black) , and other sources 
We classify as black hole candidates those X-ray 
that are resolved by Chandra as a single point, and 
Type A PDS at a 0.3-10 keV luminosity > 4x 10 37 
in at least one observation. The luminosity L c = 
for 1 4. Mq neutron star is indicated by a vertical, 
line. 



38.5 

A PDS 

(grey), 
sources 
exhibit 
erg s _1 

1 ^Edd 

dashed 



^80% were consistent with (3 ~0.5-l. Hence, broadly 
speaking, (3 ~1, consistent with Galactic LMXBs. 
However, our observed values of (3 are systematically 
lower than in the Galactic LMXBs; this is likely to be 
due to the low count rate. 

We present 0.3-10 keV luminosity functions (LFs) 
for the Type A PDS (in grey) in Fig. 3. In the solid LF 
several sources are represented more than once, since 
they exhibit Type A variability in more than one ob- 
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Table 5. Fits to the luminosity functions presented in 
Figs. 3-6. Column 1 gives the LF to be fitted, and the 
model used: power law (PL) or broken power law (BPL). 
Column 2 gives the relevant figure number, a and (5 are the 
two spectral indices below and above the break luminosity 
Lb. Prob. is the probability that the model is a good fit to 
the LF. Numbers in parentheses are 90% confidence level 
uncertainties on the last digit. 



LF a 


Fig 


. Model 


a 


P 


L b *> 


Prob. c 


A Tot 


3 


PL 


0.52(8) 






0.01 


A Tot 




BPL 


0.34(9) 


1.8(5) 


4.4(7) 


0.81 


A BH 


4 


PL 


0.40(10) 






2xl0~ 4 


A BH 




BPL 


0.05(12) 


1.5(5) 


4.6(7) 


0.97 


A NC 


4 


PL 


0.8(2) 






0.07 


A NC 




BPL 


0.5(2) 


3.6(19) 


2.1(4) 


0.89 


B Tot 


3 


PL 


0.55(2) 






0.09 


B Tot 




BPL 


0.30(2) 


0.97(7) 


3.6(3) 


0.68 


B BH 


6 


PL 


1.7(11) 






0.93 


B BH 




BPL 


2.8(2) 


1.9(8) 


4.3(17) 


0.52 


B NC 


6 


PO 


0.53(7) 






0.09 


B NC 




BPL 


0.36(10) 


0.8(2) 


5(2) 


0.56 



a Type A (A) and Type B/flat (B) LFs were modeled for 
the total (Tot), black hole (BH) and non-classified (NC) 
populations. 

b 0.3-10 keV luminosity normalised to 10 37 erg s _1 
c The goodness of fit was obtained using the method out- 
lined in (Crawford et al. 1970) 

servation; in the dotted LF each source is represented 
only once, by its highest luminosity Type A variabil- 
ity. The dashed vertical line indicates a luminosity of 
4xl0 37 erg s _1 , which we consider to be the maxi- 
mum luminosity for Type A variability in a neutron 
star (Sect. 5). We classify as a black hole candidate 
any source that exhibits Type A variability at a 0.3- 
10 keV luminosity >4xl0 37 erg s -1 in at least one 
observation. The slope of the LF is considerably shal- 
lower below 4x10 erg s _1 than above, which may 
indicate two populations, neutron star Type A's at 
low luminosity and black hole Type A's at high lumi- 
nosity. 

We modelled the Type A LF with two models: a 
power law and a broken power law. We ignored the 
Type A variability of S56, as this appears to be a 
supersoft source (See Sect. 7.5). The fitting was per- 
formed with the Sherpa modelling package, using 
Cash statistics (Cash 1979). We obtained the goodness 
of fit using the method outlined in Crawford, Jauncey, 
& Murdoch (1970) for unbinned data without uncer- 
tainties. 

We present the results of fitting the LF in Table 5. 
Uncertainties are quoted at the 90% confidence levels. 
We find that the Type A LF cannot be described by 
a single power law; the best fit has a 1% probability 
of successfully describing the data. However, it is well 
described by a broken power law where the spectral 
index changes from 0.34±0.09 to 1.8±0.5 as the lu- 
minosity exceeds 4.4±0.7xl0 37 erg s _1 in the 0.3-10 



keV band. This break may indicate the boundary be- 
tween neutron star and black hole populations, lending 
support to our classification of sources as black hole 
LMXBs if they exhibit Type A variability at 0.3-10 
keV luminosities > 4xl0 37 erg s _1 . 

In determining our observed luminosities, we as- 
sume that the emission is isotropic. Instead, the emis- 
sion could be anisotropic and beamed in our direction. 
In this case, the range of high luminosity Type A vari- 
ability could be due to varying beaming factors rather 
than variation in primary masses. This possibility is 
further explored in Sect. 7.4.1. 

In Fig. 4, we separate the Type A LF into black 
hole candidates and other X-ray sources, most likely 
neutron star LMXBs. However, this non-classified 
population may contain hidden black holes where I 
<C l c . Here, the dashed line represents 2xl0 37 erg s _1 , 
which is our value of L c for a 1.4 M Q neutron star. 
Fits to the black hole and non-classified Type A LFs 
are presented in Table 5. 

The Type A LF for the black hole candidate pop- 
ulation appears well described by a simple power law 
except at the low luminosity end where we are neces- 
sarily incomplete. However, a power law fit is rejected 
because some of the black hole candidates also ex- 
hibited Type A variability at luminosities <4xl0 37 
erg s _1 . Instead, the LF is well described by a bro- 
ken power law with the spectral index changing from 
0.05±0.12 to 1.5±0.5 at a luminosity of 4.6±0.7x 10 37 
erg s _1 . 

The Type A LF for the non-classified population is 
best described by a broken power law with the spectral 
index changing from 0.5±0.2 to 3.6±1.9 at a luminos- 
ity of 2.1±0.4xl0 37 erg s _1 . Interestingly, this break 
agrees well with our predicted maximum luminosity 
for Type A variability in a 1.4 M Q neutron star, sug- 
gesting that this non-classified population is indeed 
dominated by neutron star LMXBs. However, the best 
fit power law model has a goodness of fit probability 
of 7%, and hence cannot be rejected. 

We also note that ~45% of the unidentified Type A 
PDS were observed at less than 1.0 x 10 37 erg s _1 ; this 
strongly argues against the presence of a bias towards 
higher luminosity (i.e. signal to noise) that would oc- 
cur if the Type A phenomenon was due to statistical 
rather than intrinsic noise in the systems. 

7.1.2. Type B / Flat PDS 

In Fig. 5, we present luminosity functions of Type 
B and flat PDS, normalised by the total number in 
each sample. KS-testing gave a probability of 47% that 
the two LFs were drawn from the same distribution. 
Hence, we consider the Type B and flat PDS to be 
consistent with being from the same distribution, sug- 
gesting that most of the flat PDS in our sample in fact 
come from high accretion rate LMXBs. 

The combined luminosity function of flat and Type 
B PDS is also shown in Fig. 5. This total LF for Type 
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B and flat PDS is compared with the Type A LF 
in Fig. 3. Best fit power law and broken power law 
models for the total Type B/flat LF are presented in 
Table 5. The best fit power law model provides an ac- 
ceptable fit (9% probability of a good fit); however, 
a broken power law is prefered (68% probability of 
a good fit). The spectral index of the best fit broken 
power law model changes from 0.30±0.02 to 0.97±0.07 
at 3.6±0.3xl0 37 erg s^ 1 . We see from Fig. 3 and 
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Fig. 5. Normalised, fractional luminosity functions for 
Type B (light grey) and flat (dark grey) PDS. KS test- 
ing shows that these two populations have a 47% chance 
of being drawn from the same distribution. The similar- 
ity between these supports the idea that the flat PDS are 
consistent with Type B PDS. The combined fractional lu- 
minosity function (black) of Type B and flat PDS is also 
given. 
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Fig. 6. Cumulative luminosity functions for Type B/flat 
PDS in black hole candidates (black), and other sources 
(grey). The axes are log-scaled. The lowest luminosity 
Type B/flat PDS for the black hole candidates is consis- 
tent with ~4xl0 37 erg s" 1 (in S63, Table A. 2). This result 
supports l c ~0.1. We note that there are no observations 
of Type B variability in the black hole candidates at 0.3- 
10 keV luminosities >10 erg s _1 . This is not unexpected, 
since Galactic black hole LMXBs are observed in the steep 
power law state at high luminosities, and exhibit Type A 
variability. 



Table 5 that the Type A and Type B/flat LFs are dis- 
tinct. Hence we can be confident that the break in the 
Type A LF is not some artifact of the XMM-Newton 
observations. 

Figure 6 shows the combined Type B/flat luminos- 
ity functions for the black hole candidate and uniden- 
tified populations. Even though the black hole candi- 
dates were selected for their Type A PDS, no Type 
B/flat PDS is observed below ~4xl0 37 erg s _1 in this 
black hole population. Hence these results support l c 
-0.1. We see from Table 5 that the Type B/flat LF for 
the black hole candidates is best described by a power 
law with spectral index 1.7± 1.1. The Type B/flat 
LF for the non-classified population is consistent with 
that of the total population. 

In Fig. 6 we see four observations of Type B/flat 
PDS in the unidentified population at 0.3-10 keV lu- 
minosities in excess of 3 x 10 38 erg s _1 , well above LEdd 
for a 1.4 M Q neutron star. Three of these are observa- 
tions of the Z-source S26; however, the emission is lo- 
cally sub-Eddington, with the neutron star contribut- 
ing only 7-15xl0 37 erg s _1 (Barnard et al. 2003c). 
The remaining emission was contributed by the ADC, 
which is probably extended (see Sect. 2), and hence 
locally sub-Eddington also. The other high luminosity 
observation was from the supersoft source S35, which 
may be beamed (see Sect. 7.5). 

7.1.3. Variability of the PDS 

The PDS of Galactic LMXBs are characterised by 
their r.m.s. variability, as detailed in Sect. 4. We classi- 
fied our M31 X-ray sources using the fractional r.m.s. 
variability amplitude, F vav , accounting for the uncer- 
tainties in the measurements (see e.g. Vaughan et al. 
2003). We then compared the variability of Type A 
and Type B/flat PDS in our sample with the known 
r.m.s variability of Galactic LMXBs. 

The r.m.s. variability of our sample is measured 
in the 0.001-0.1 Hz range, while r.m.s. variability 
of Galactic LMXBs is often measured in the 0.1-10 
Hz range (see e.g. McClintock & Remillard 2004). 
The PDS and r.m.s. variability properties of Galactic 
LMXBs are considerably better known than those of 
our sample; hence, we estimated the r.m.s. variabil- 
ities of Galactic LMXBs in the 0.001-0.1 Hz range. 
The r.m.s. variability of a lightcurve is related to the 
integrated PDS (see e.g. Vaughan et al. 2003). Hence, 
if V\ Q and Vhi represent the source variability in the 
0.001-0.1 Hz and 0.1-10 Hz ranges respectively, and 
A\ Q and A^i represent the area under the PDS in the 
same frequency ranges, then 

— ^ (21 
V hi ~ A hi ' {) 

To estimate Vi /Vhi, we created artificial Type A 
PDS with the spectral index 7 changing from to 
1 at break frequencies of 0.01 and 0.07 Hz, span- 
ning the frequency range of observed break frequen- 
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Fig. 7. Fractional r.m.s. variability, F var vs. 0.3-10 keV 
luminosity for Type A (top panel) and Type B/flat (bot- 
tom panel) PDS. The axes are log-scaled. Shaded regions 
show estimated r.m.s. variabilities for Galactic Type A and 
Type B PDS in the 0.001-0.1 Hz range. A power law fit to 
the Type A data is shown in both panels, for comparison 
of the variabilities of Type A and Type B/flat PDS. 



cies within 90% confidence limits (see Table A. 2). We 
found Vio/Vhi to be —0.7 and 0.3 respectively. This 
suggests that we are most sensitive to systems with 
low break frequencies. We also created Type B PDS, 
using a simple power law model with 7 = 1.0 and 
1.5, spanning the range of spectral indices for Type 
B variability in Galactic LMXBs, giving V\ /Vhi —1.1 
and ~13 respectively. However, we will need to exam- 
ine the 0.001-0.1 Hz PDS of Galactic LMXBs in detail 
at a later date. 

We present plots of F var vs. 0.3-10 keV luminosity 
for the classified PDS in Fig. 7; the results for Type 
A PDS are shown in the top panel, and Type B/flat 
results are shown in the bottom panel. The error bars 
represent the 90% confidence limits. In each panel, the 
shaded area shows the range of r.m.s. variability for 
Galactic LMXBs in the 0.001-0.1 Hz range, estimated 
as above. We also plot in each panel the best fit power 
law model to the Type A data. 

We see that the variabilities of most of the PDS 
are consistent with their Galactic counterparts. The 
Type B PDS are systematically more variable than 
Galactic counterparts with 7 = 1.0; however, they are 
consistent with the known range of 7 for Galactic Type 
B variability: 1.0 < 7 < 1.5 (van der Klis 1994). The 
Type A variability is broadly characterised by F V&T 
-0.42 L^ 059 , where L 37 is the 0.3-10 keV luminosity 
in units of 10 37 erg s _1 . We see from the bottom panel 
of Fig. 7 that F var <0.42 L^?' 59 for most Type B PDS; 
hence our Type A PDS tend to be more variable than 
our Type B PDS as expected, but this result is not 
strongly significant. 



7.2. The bright globular cluster X-ray sources 

There are thirteen bright X-ray sources that have been 
identified in Galactic globular clusters. Twelve of these 
have been identified as neutron star LMXBs, while 
the thirteenth has not been classified (in't Zand et al. 
2004, and references therein). Hence the 14 globular 
cluster X-ray sources in our sample are expected to 
be LMXBs containing —1.4 M Q neutron stars. This 
makes their observed spectral and timing behaviour 
of particular interest. 

We summarise the spectral and timing properties 
of the globular cluster sources in Table 6. We stress 
that S47 is resolved into two unrelated sources by 
Chandra; the southern source is associated with the 
globular cluster, while the northern source is transient 
and contains a black hole candidate (Barnard et al. 
2004). Our results for the globular cluster sources 
are consistent with Type A variability being observed 
at lower luminosities than Type B. Furthermore, our 
SEDs are well described by power laws with photon 
index —1.5-2.1 when Type A variability is observed, 
consistent with LMXBs in the low/hard state as clas- 
sified by McClintock & Remillard (2004). In this re- 
spect, the globular cluster X-ray sources behave in a 
way consistent with Galactic neutron star LMXBs. 

In Fig. 8 we show PDS type vs. luminosity for ob- 
servations of the globular cluster sources where the 
PDS is classified; Sll is omitted because it showed 
no classified variability. As before, we assume that l c 
—0.1 in the 0.3-10 keV band. A dashed vertical line 
represents a luminosity of 2xl0 37 erg s _1 , i.e. —10% 
Eddington for a 1.4 M neutron star. Out of the 12 
globular cluster X-ray sources in Fig. 8, 3 systems 
show more than one type of variability (S4, S21, and 
S43), while 6 systems are observed at luminosities 
near L c -2xl0 37 erg s" 1 (S4, S15, S21, S43, S48, 
S49). Hence we can only estimate Z c with these data. 
Ten of the twelve sources behaved as expected for 
LMXBs with a 1.4 M Q neutron star primary. However, 
S4 and S43 appear at face value to be discrepant. We 
discuss how these sources may fit into our scheme be- 
low. 

Most striking are the two observations of Type A 
variability in S43, at luminosities >5xl0 37 erg s~ 1 . 
A third observation of S43 has a flat PDS, at a lu- 
minosity that is consistent within errors with being 
greater than the Type A luminosity in S43. Hence, 
S43 is consistent with l c —0.1 if it contains a black 
hole rather than a neutron star. Although no Galactic 
black hole LMXBs have been identified in globular 
clusters, we note that Angelini et al. (2001) have re- 
ported a possible globular cluster black hole binary 
in the elliptical galaxy NGC 1399, with a 0.3-10 keV 
luminosity of 5xl0 39 erg s _1 . 

Other possibilities are explored by Di Stefano et al. 
(2002) when discussing an M31 globular cluster X- 
ray source in the disc with a 0.2-7 keV luminos- 
ity of 2-6 xlO 38 erg s _1 . Firstly, the globular cluster 
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Table 6. Luminosities, SED shapes and PDS for each of the 14 bright globular cluster X-ray sources in Observations 
1-4. Ln is the 0.3-10 keV luminosity in Observation TV, in units of 10 37 erg s _1 ; Sn is the photon index for the best 
to the SED of a power law model, and Tjv is the PDS type. ". . ." indicates that no good SED was obtained; either the 
source was not observed with the PN, or the SED had negative counts. 
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a S47 is resolved into two sources by Chandra, the globular cluster X-ray source to the south, and a black hole 
candidate to the north. 

b The spectra of S55 in observations 1 and 4 are particularly hard because non-dip and dip spectra are mixed. 



could contain two or more X-ray sources. White & 
Angelini (2001) discovered one globular cluster with 
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Fig. 8. PDS type and 0.3-10 keV luminosity for observa- 
tions of globular cluster X-ray sources where the PDS type 
is classified. A dashed line represents 2xl0 37 erg s _1 , i.e. 
10% Eddington for the assumed 1.4 Mq primary. The y- 
axis is used to separate the results of the different sources 
and has no physical meaning. Circles, squares and trian- 
gles represent observations of Type A, Type B and flat 
PDS respectively. We expect Type A PDS at luminosities 
below, and Type B or flat PDS above, 2xl0 37 erg s~\ 
S43 is inconsistent with L c — 2xl0 37 erg s _1 , as it ex- 
hibited Type A variability at a 0.3-10 keV luminosity of 
5-6xl0 37 erg s~\ S4 exhibited Type B and flat PDS at 
lower luminosities than expected; this may be due to a low 
mass primary, or hysteresis. 



a second bright X-ray source out of a population 
of 13 Galactic X-ray bright globular clusters, hence 
it is reasonable to expect that ~1 of the 14 X-ray 
bright globular clusters in our sample contains mul- 
tiple LMXBs. Alternatively, S43 could contain a 1.4 
M Q neutron star primary, accreting ^100% helium, 
rather than hydrogen; this would raise the Eddington 
limit by a factor of 2. Finally, the observed luminosity 
could be artificially enhanced by relativistic beaming, 
or overestimated if anisotropic; this will be discussed 
in Sect. 7.4.1. 

S4 exhibited flat and Type B PDS in observations 
1 and 2 respectively, at 0.3-10 keV luminosities of 
1.3±0.3 and 1.0±0.3 xlO 37 erg s _1 . According to our 
model, Type A PDS are expected at these luminosi- 
ties. S4 is consistent with l c = 0.07 LEdd for a 1.4 
Mq neutron star primary. However, S4 would be also 
consistent with l c = 0.1 for a LMXB with a ~1.0 M Q 
neutron star primary. There is already observational 
evidence for systems with low mass neutron stars. An 
example is the eclipsing binary pulsar SMC X-l, which 
is discussed next. 

Van der Meer et al. (2005) measured the mass of 
the primary in SMC X-l to be 1.05±0.09 M , using 
high resolution radial velocity measurements, and as- 
suming that the secondary filled its Roche lobe. Val 
Baker, Norton, & Quaintrell (2005) also measured the 
mass of the neutron star in SMC X-l. They obtained 
a maximum neutron star mass by assuming a Roche 
lobe filling secondary, and also obtained a minimum 
mass by assuming an inclination of 90°. In each case 
they corrected for X-ray heating of the secondary, un- 
like van der Meer et al. (2005). Val Baker et al. (2005) 
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Table 7. Luminosities, SED shapes and PDS in Observations 1-4 for each of the 14 black hole candidates, either 
exhibiting Type A variability at 0.3-10 keV luminosities > 4 x 10 37 erg s— 1 in any observation, or identified by some 
other method. Ln is the 0.3-10 keV luminosity in Observation N, in units of 10 37 erg s _1 ; Sn is the photon index for 
the best to the SED of a power law model, and Tn is the PDS type. ". . ." indicates that no good SED was obtained; 
either the source was not observed with the PN, or the SED had negative counts. 
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a Previously classified as Type A by Barnard et al. (2003b), using different good times, frequency ranges and PDS 
classifications. 

^Located within a globular cluster 

c Chandra resolved S47 into two sources; the information presented here is based on combined data from Chandra and 
XMM-Newton observations (Barnard et al. 2004) 
identified as a possible black hole by Williams et al. (2004a) 

e The high spectral indices presented in this table do not indicate the best fit models. 
f S63 was not in the EPIC field of view for Observations 1 and 3. 



obtained upper and lower mass limits for the neutron 
star of 1.01±0.10 and 0.73±0.08 M respectively be- 
fore correction; this upper limit agrees well with the 
mass found by van der Meer et al. (2005). After cor- 
rection, maximum and minimum masses of 1.21±0.10 
and 0.91±0.08 M Q were obtained. Hence, the low lu- 
minosity Type B and flat PDS observed in S4 may be 
due to a low mass primary. 

Alternatively, it is possible that S4 is hysteretic, 
with the transition from Type B to Type A occurring 
at a lower luminosity than the transition from Type A 
to Type B, like the Galactic neutron star LMXB Aql 
X-l (see e.g. Maccarone & Coppi 2003). 

7.3. X-ray sources associated with supernova 
remnants 

Two sources, S38 and S58, are associated with SNRs, 
but also exhibit Type A PDS. In fact, the SNR of S58 
has been resolved by Chandra (Kong et al. 2002b). 
Only a handful of accreting compact objects have 
been associated with SNRs: SS 433 (e.g. Dubner et al. 
1998), MF16 in the spiral galaxy NGC 6946 (e.g. 
Roberts & Colbert 2003), XTE J0111. 2-7317 in the 
SMC(Coe et al. 2000) and RX J050736-6847.8 in the 
LMC (Chu et al. 2000). None of these are thought to 
be LMXBs. 

Monte-Carlo simulations of LMXB formation have 
shown that neutron star LMXBs have an average sys- 



temic velocity of 180±80 km s _1 (sec e.g. Brandt & 
Podsiadlowski 1995). The X-ray remnant of S58 is 40 
pc in diameter, and Kong et al. (2002b) estimate the 
age of the remnant to be ~3, 000-20, 000 yr. Hence, the 
neutron star is unlikely to have travelled further than 
~10 pc. The extraction region for S58 has an angu- 
lar size of 40 , corresponding to ~150 pc; hence any 
LMXB formed in the supernova would be within our 
extraction region. We further discuss S58 in a separate 
paper (Williams et al. 2005). 

7.4. The black hole candidates 

We found 12 sources that exhibited Type A variability 
at 0.3-10 keV luminosities >4x 10 37 erg s _1 in at least 
one observation. In addition, the PDS of S14 were pre- 
viously identified as Type A in Observations 1, 3 and 
4 by Barnard et al. (2003b), using different good time 
intervals and frequency ranges; also S51 was identified 
as a black hole candidate by Williams et al. (2004a) 
from its transient behaviour. All these sources are seen 
as a single point source by Chandra. 

The 14 systems are listed in Table 7; the luminos- 
ity, spectral shape and PDS type are given for each of 
the four XMM-Newton observations of each system. 
The spectral shape is illustrated by the photon index 
of the best fit power law model. 

We have already published results for two of our 
black hole candidates. S14 is the prototype for our 
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new method for identifying black hole candidates, and 
is discussed in Barnard et al. (2003b). S47 is a point 
source in the XMM-Newton observations, but was re- 
solved by Chandra into two sources; we identified the 
northern source as a black hole candidate after com- 
bining data from these XMM-Newton observations 
with 35 Chandra observations (Barnard et al. 2004). 

7.4.1. Could the high luminosity Type A variability 
be due to beaming? 

Even if we assume that l c ~0.1, systems that exhibit 
Type A variability at high luminosities could contain 
a neutron star, or even a white dwarf, if the emis- 
sion were highly beamed or anisotropic. Such beam- 
ing could arise from the relativistic boosting of jet 
emission, or from reprocessing of X-rays by a warped 
inner disc (c.f. Pringle 1996; Foulkes et al. 2006). 
Alternatively, the disc could simply have a lower op- 
tical depth to scattering in some directions(e.g. along 
the rotational axis) than in others (see e.g. King et al. 
2001). 

Narayan & McClintock (2005) studied 10 Galactic 
black hole candidates to determine whether their qui- 
escent luminosities were dependent on their inclina- 
tions; i.e., whether the quiescent emission was signif- 
icantly beamed. However, they found no significant 
trend to suggest that relativistic beaming was signif- 
icant. They compared the observed luminosities and 
inclination angles with beaming models for different 
Lorentz factors, 7, and found that 7 <1.24 at the 90% 
confidence limit, if the X-ray emission originates from 
a jet. Similarly, Maccarone (2003) found that 7 ^1.4 
in the X-ray emitting gas in GRO J1655— 40 during 
the low/hard state. With such low Lorentz factors ob- 
served, we therefore conclude that the high luminosity 
Type A variability is unlikely to be due to relativistic 
Doppler boosting. 

We cannot rule out the possibility that any of our 
sources are emitting anisotropically. Assuming l c ~0.1, 
and a 1.4 M Q primary, we define a hypothetical beam- 
ing factor for a source, bh, as the highest luminos- 
ity for Type A variability in the source, divided by 
2xl0 37 erg s _1 . If the anisotropic emission were asso- 
ciated with the disc, then it would be feasible for the 
beam to precess. Hence, a source that varies greatly 
between observations may be emitting in precessing, 
anisotropic beam. S56 is one such source, and is fur- 
ther discussed in Sect. 7.5. 

7.4.2. X-ray properties of the black hole candidates 

Thirteen out of the fourteen sources listed in 
Table 7 have SEDs consistent with Galactic black hole 
LMXBs. However, S56 is a supersoft source, with a 
temperature of ~20-30 eV, and cannot be considered 
as a candidate black hole LMXB in the low/hard state. 
We discuss the possible nature of S56 in Sect. 7.5. We 



now briefly describe the properties of each black hole 
candidate in turn. 

• S7 We observed Type A variability from S7 
at a 0.3-10 kcV luminosity of 8±2xl0 37 erg s" 1 in 
Observation 4. However, S7 exhibited Type B vari- 
ability at 3.7±0.3 and 5.3±0.3 xlO 37 erg s" 1 in obser- 
vations 1 and 3 respectively. The power law compo- 
nent of the SED of S7 in Observation 4 has a spectral 
index of 2.47±0.12; hence, S7 appears to be in the 
steep power law state, as defined by McClintock & 
Remillard (2004). 

Assuming l c ~0.1, then the low luminosity Type 
B variability observed in Observation 1 indicates that 
the primary is relatively low mass. However, since S7 
was not seen in either the Einstein or ROSAT surveys, 
it could be a recurring transient, and in decay from 
outburst during Observations 1 and 3. In this case, 
the transition from Type B to Type A during the 
decay of the outburst would occur at a significantly 
lower luminosity than L c . This scenario would, how- 
ever, require at least two outbursts in the 18 months 
between Observation 1 and Observation 4. 

If the emission from S7 were beamed, then the 
beaming factor bh would be <;3 if the primary is 1.4 
M Q neutron star. The brightest observation of S7 is 
only ~50% brighter than the faintest, so we don't con- 
sider S7 to be a strong candidate for beaming. 

• S12 Observation 4 revealed Type A vari- 
ability in S12 at a luminosity of 6.0±0.3 xlO 37 erg 
s _1 . The other observations of S12 had luminosities 
4.6-5.0 xlO 37 erg s _1 , but no PDS were classified. 
However, all the SEDs from S12 were consistent with 
low/hard state SEDs (i.e., consistent with T ~1.5- 
2.1). 

• S14 There were three differences in the anal- 
ysis of S14 between the current paper and Barnard 
et al. (2003b). Firstly, we were able to use the whole 
of Observation 1, rather than a shorter good time in- 
terval, as the background flaring did not significantly 
contribute to the 0.3-10 keV flux. Secondly, our PDS 
were constructed from 512 bins of 2.6 s. Finally, and 
most importantly, we allowed a, the spectral index of 
the PDS below the break frequency, to vary. In ob- 
servations 1, 3 and 4, a was found to be 0.11±0.01, 
0.24±0.03 and 0.27±0.04 respectively; hence the PDS 
were not well fitted by models where a = 0. In addi- 
tion, the breaks occurred at high frequencies, so that 
the PDS were difficult to classify using frequencies 
<0.1 Hz. 

We showed in Barnard et al. (2003b) that a two 
component model was required to fit the Observation 
4 SED of S14, comprising a 1.2-2.0 keV black- 
body/disk blackbody component and a power law 
component with T — 2.4-3.4. We argued incorrectly 
that the SEDs of Observations 1 and 3 must also have 
two components; while the SEDs from Observations 
3 and 4 are very similar in shape and flux, the SED 
of Observation 1 is well described by a single power 
law with T <~1.7. We now conclude that S14 was in 
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the low/hard state in Observation 1 and the very high 
state in Observations 3 and 4. 

• S19 Type A variability was observed at 
9.7±0.5x 10 37 erg s _1 in Observation 4 , while Type B 
variability was observed at 3.3±0.3 and 5.1±0.5 x 10 37 
erg s^ 1 in Observations 1 and 3 respectively. The SED 
for Observation 4 is fitted by T = 2.47±0.09, which is 
consistent with the very high (steep power law) state. 
However, the low Type B luminosity in Observation 1 
indicates a low primary mass, perhaps —3 M Q if l c 
-0.1. 

• S37 We see Type A variability from S37 in 
Observations 3 and 4, at luminosities of 8.6±0.3 and 
8.7±0.3 xlO 37 erg s" 1 respectively. The SEDs of these 
observations fit with T —1.8-1.9, consistent with the 
low/hard state. A flat PDS is observed in Observation 
2, at a luminosity of 8.9±0.7xl0 37 erg s -1 , while the 
PDS from Observation 1 is unclassified. S37 is there- 
fore consistent with being a canonical black hole. 

• S40 Type A variability is seen in Observation 3 
with a SED fitted by T = 1.88±0.09, at a luminosity 
of 4.9±0.3xl0 37 erg s" 1 . Additionally, a flat PDS was 
observed in Observation 4 at 4.3±0.2xl0 37 erg s _1 . 
S40 is therefore consistent with having a black hole 
primary, with a relatively low mass. 

• S43 Observations 2 and 3 revealed Type A 
variability with SEDs consistent with the low/hard 
state at luminosities of 5.1±0.6 and 5.5±0.3 xlO 37 
erg s _1 respectively. In addition, S43 had a flat PDS 
in Observation 1, at a luminosity of 4.8±0.4x 10 37 erg 
s _1 . However, S43 is associated with a globular cluster 
(see Sect. 7.2). The maximum variation in luminosity 
—20%, while bh ^3; it may be beamed, but we have 
no evidence to show it. 

• S47 Chandra resolves S47 into two sources. The 
northern source is transient, and contains the black 
hole candidate, while the southern source is associated 
with a globular cluster, (see Barnard et al. 2004). 

• S51 Williams et al. (2004a) identified S51 as a 
black hole candidate, from its transient behaviours ob- 
served in Chandra and HST observations. Our XMM- 
Newton observations reveal little extra information 
about the system. S51 exhibits Type A variability 
in Observation 3, hence it is probably disc accreting. 
However, the 0.3-10 keV luminosity derived from our 
best fit to the SED is only 4±4xl0 36 erg s _1 , so we 
cannot use high luminosity Type A variability to clas- 
sify it is a black hole candidate. 

• S53 In Observations 1-3, S53 is barely detected, 
and the PDS are unclassified. In Observation 4, a Type 
A PDS is observed at 4.2±0.2xl0 37 erg s _1 , while the 
SED is fitted with V = 1.83±0.08, consistent with the 
low hard state. We therefore consider S53 to be a M31 
counterpart to the Galactic black hole transients. 

• S59 Type A variability was seen in Observations 
2 and 3, at luminosities of 9.1±1.5 and 6.8±0.7 xlO 37 
erg s^ 1 respectively. A flat PDS was observed at 
7.8±1.0xl0 37 erg s _1 in Observation 1, consistent 
with our black hole model. However, the SEDs are 



fitted with T —0.9-1.1; i.e. the SEDs are too hard for 
the low/hard state in LMXBs. Instead, S59 could be 
a black hole HMXB. 

• S60 The SEDs of S60 are also too hard for 
LMXBs in the low/hard state: T -0.8-1.1. Type A 
variability was observed at luminosities of 6. 8 ±1.2 and 
5.8±0.3 xlO 37 erg s -1 ; no other PDS were classified. 
S60 may also be a HMXB. 

• S63 The source was only in the field of view in 
Observations 2 and 4. Type A variability was seen in 
both, at luminosities of 3.1±0.6 and 6.4±0.5 xlO 37 
erg s -1 respectively. However, the SEDs of S63 in 
Observations 2 and 4 are too soft for low/hard state 
spectra, suggesting that we are observing S63 in the 
steep power law state. The steep power law state is 
sometimes observed at the transition from Type A to 
Type B; hence, one of these luminosities may corre- 
spond to L c in the 0.3-10 keV band. If S63 is hys- 
teretic, like the Galactic black hole LMXBs, then L c 
could be 6.4±0.5xl0 37 erg s _1 . Otherwise, L c could 
be as low as 3.3±0.8xl0 37 erg s -1 ; this lower L c would 
indicate a primary mass of —3 M Q . 

7.4.3. A new population of persistently bright 
black hole LMXBs? 

Perhaps the most surprising result of this survey 
is that only 3 out of the 13 black hole candi- 
dates are transient, and 8 of the remaining 10 are 
most likely LMXBs. By contrast, all of the confirmed 
Galactic black hole LMXBs are transient. Indeed S14, 
S19, S40, S59 and S63 are present in the 1979-1980 
Einstein surveys (Trinchicri & Fabbiano 1991), the 
1991 and 1992 ROSAT surveys (S97, S01), as well as 
all pointed Chandra ACIS (W04) and XMM-Newton 
observations. The SED of S59 suggests it may be a 
HMXB, but the other four sources appear to be black 
hole candidates in persistently bright LMXBs. 

Alternatively, these sources could instead be ex- 
hibiting long duration, —20 year, outbursts, similar to 
Galactic black hole LMXBs GX 339-4 (e.g. Zdziarski 
et al. 2004) or GRS 1915+105 (e.g. Belloni et al. 2000). 
However, the first extra-solar X-ray source was discov- 
ered only —40 years ago (Giacconi et al. 1962); hence 
one might well ask if any LMXB is truly persistent. 

The existence of persistently bright LMXBs is in 
fact predicted by theory. Evolutionary theory pre- 
dicts a substantial population of short period black 
hole LMXBs with unevolved companions that would 
exceed the critical mass transfer rate for stable disc ac- 
cretion (see King et al. 1997, and references within). 
Hence large numbers of persistently bright black hole 
LMXBs would be expected. 

Disc accretion at low mass transfer rates is un- 
stable; however, in LMXBs, disc irradiation plays a 
role (van Paradijs 1996; King et al. 1996). The critical 
accretion rate that stabilises the disc, Merit' * s propor- 
tional to where Ru is the disc radius; hence, M^it 
oc P 4 / 3 by Kepler's laws, where P is the orbital period. 
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The fact that Galactic black hole LMXBs are transient 
where neutron star LMXBs with a similar period are 
persistent suggests that irradiation is weaker in black 
hole LMXBs (King et al. 1997). 

Using the assumptions given in King et al. (1997), 

M£ t ~ 2.86 x lO-^V^ffXyr- 1 (3) 

where m\ is the primary mass and rri2 is the secondary 
mass in solar units, and Ph is the orbital period in 
hours (see also Dubus et al. 1999). King et al. (1997) 
calculated M over Ph ~l-8 hr for black holes with 
masses 2, 5 and 10 M Q , for comparison with M"\ t . 
Both the 2 and 5 M Q black holes were sufficiently ir- 
radiated for a stable disc in the range ~3-8 hr. For the 
10 M Q black hole, M < M' c " t for all periods. Hence, 
LMXBs with lower mass primaries are more likely to 
be persistent. However, due to uncertainties in their 
assumptions, King et al. (1997) acknowledge that even 
10 Mq black holes could be persistent in this range 
of periods, while long-period LMXBs are always ex- 
pected to be transient (King 2000). 

7.4.4. Multi-wavelength observations of the black 
hole candidates 

Of the thirteen black hole candidates, only S43 and 
S53 have been observed at other wavelengths. S43 is 
associated with the globular cluster Bo 144 (Battistini 
et al. 1987), and has been observed in UV the optical 
and infrared (Galleti et al. 2004); however, the coun- 
terpart to S43 will be lost in the cluster, so we will 
discuss S43 no further. Williams et al. (2004a) report 
results from Chandra/HST observations of S53, using 
the designation "r3-16". They find a roughly persis- 
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Fig. 9. PDS of S53 from Observation 4, using 192 inter- 
vals of 64 bins. We show the best fits to the PDS using 
power law (solid) and broken power law (dashed ) models. 
For the best fit power law and broken power law models, 
X 2 /dof = 35/9 and 11/8 respectively; the broken power 
law is formally acceptable, while the power law is not. We 
also show the PDS of the background, in grey with open 
circles. It is clear that the source PDS is not dominated 
by background variability. 



tent U band counterpart to r3-16 that resembles an 
unresolved binary; the U magnitude was measured to 
be 20.82±0.06 in 2001, August and 2001, December, 
and 21.11±0.02 in 2002, January. Measurements from 
the WFPC2 instrument on board HST suggest a sepa- 
ration between the two components of 0.82±0.08 pc, 
if located in M31. However, ground based measure- 
ments using the B and V band data from the Local 
Group Survey of Massey et al. (2001) yield a sepa- 
ration of 3.9±0.6 pc. Williams et al. (2004a) suggest 
that S53 could be either a blend of foreground stars, 
a cataclysmic variable (CV) at a distance ^1 kpc, or 
a background AGN. They note that the X-ray SED is 
harder than usual for a CV, and that the SED and 
linc-of-sight absorption are more in keeping with a 
background AGN. However, the Type A PDS exhib- 
ited by S53 in Observation 4 means that it cannot be 
a background AGN, and is unlikely to be a blend of 
foreground stars. 

Furthermore, the 37W radio survey of M31 
(Walterbos et al. 1985) revealed a 1.1 mJy radio source 
(37W153) just 0.6" from the Chandra position of 
S53; the survey was conducted at a frequency of 1412 
MHz and had a 23" x 35" resolution. Radio emis- 
sion is common in black hole LMXBs in the low state 
(e.g. Fender et al. 2004), but is rare in CVs. Claims 
have been made for detecting radio emission from sev- 
eral magnetics CVs (e.g. Chanmugam 1987; Pavclin 
et al. 1994), as well as in novae, dwarf novae, and re- 
current novae (Chanmugam 1987). However, the only 
confirmed radio detection is from the peculiar mag- 
netic CV AEAquarii (e.g. Meintjes & Venter 2003). 

Cygnus X-3 is one of the brightest X-ray binary 
radio sources; it exhibits flares up to ^20 Jy (see 
e.g. Waltman et al. 1989). Predehl et al. (2000) es- 
timate the distance of Cygnus X-3 to be ~9 kpc, us- 
ing Chandra observations of its scattering halo. If we 
place Cyg X-3 in M31, the radio flares can reach ^2.5 
mJ. Hence, it is feasible for S53 to be a M31 analogue 
to Cyg X-3. 

S53 is unlikely to be a magnetic CV since its Type 
A variability (Fig. 9) indicates disc accretion. Williams 
et al. (2004a) found no clear correlation between the 
X-ray flux of S53 and the U band flux of its coun- 
terpart; hence, the association between the UV and 
X-ray sources may be coincidental. 

7.5. The supersoft sources 

There are four sources in the survey that have been 
classed as supersoft: S35, S56, S57 and S58. Many 
supersoft sources are thought to be accreting white 
dwarf binaries, with sufficiently high mass transfer 
rates for stable hydrogen burning (van den Heuvel 
et al. 1992). However, S58 is a supernova remnant, 
and discussed in Sect. 7.3. 

Blackbody radiation models assume that the radi- 
ation source is optically thick at all energies. However, 
the opacities of model white dwarf atmospheres are 
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dominated by only a few ions in the temperature range 
typical of supersoft sources; as a result, the X-ray 
spectra deviate significantly from blackbody emission 
(Hcise et al. 1994). A blackbody fit to a white dwarf 
atmosphere will significantly over-estimate low-energy 
flux and radius (Heise et al. 1994). Indeed, blackbody 
fits to the spectra of S35 and S56 could be formally 
rejected in Observations 3 and 4, with best fit x 2 /dof 
>2. 

Detailed modelling of the supersoft sources with 
white dwarf atmosphere models is beyond the scope of 
this paper. However, neutron star atmosphere (NSA) 
models give qualitatively similar SEDs, so we applied 
NSA models to the supersoft sources as a first approx- 
imation. 

S35, S56 and S57 have broadly similar character- 
istics; all are well described by a NSA model, with 
temperatures of 20-30 eV, and apparent 0.3-10 keV 
luminosities >10 38 erg s _1 in at least one observation. 
However, S57 is a transient that exhibited —850 s pul- 
sations in Observation 1 (Osborne et al. 2001), while 
the 0.3-10 keV luminosity of S56 apparently varied 
over 1.5-40 xlO 37 erg s _1 , exhibiting the highest lu- 
minosity Type A variability in the survey. 

The pulsations exhibited by S57 shows that the 
emission is anisotropic. Furthermore, S56 exhibited 
Type A variability at — 4xl0 38 erg s _1 , suggesting 
&h ^20; S56 also varied by a factor of —20 between 
observations, indicating that the emission is likely 
to be anisotropic. Hence we could expect the emis- 
sion from S35 to be anisotropic also. Indeed, col- 
limated jets have been observed from the Galactic 
supersoft sources RX J0925. 7-4578 (Motch 1998), 
RXJ0513.9-6951 (Crampton et al. 1996; Southwell 
et al. 1996) and RXJ0019. 8+2156 (Tomov et al. 1998; 
Becker et al. 1998). It is therefore quite possible that 
S35, S56 and S57 are all disc accreting white dwarf 
binaries that exhibit some sort of jet. 

8. Summary 

We have conducted a survey of PDS and SEDs ob- 
tained from 4 XMM-Newton observations of the cen- 
tral region of M31. The sample consisted of the 63 
X-ray sources that exceeded a 0.3-10 keV EPIC-pn 
intensity of 0.02 count s _1 in at least one observa- 
tion. Four XMM-Newton sources were identified as 
multiple sources by Chandra: S22, S31, S33 and S47; 
in Barnard et al. (2004), we show that the northern 
source in S47 is a possible black hole. A further four 
sources are supersoft sources. The known foreground 
star, S6, showed neither Type A nor Type B PDS, 
as expected. The remaining 54 sources are consistent 
with being LMXBs; we identify 13 as black hole can- 
didates. 

Variability over a 100 s time scale is observed in 
55 out of the 63 sources. Of these, 40 sources show 
Type A variability, meaning that they cannot be back- 
ground AGN and are most likely to be LMXBs; a fur- 



ther 6 show Type B but not Type A PDS, making 
them likely accreting binaries. Hence, we identify a to- 
tal of 46 likely X-ray binaries out of 63 X-ray sources, 
from their PDS and SEDs. We find that the PDS and 
SEDs correspond very well to the spectral and timing 
states of Galactic X-ray binaries; 2 sources appear to 
be HMXBs while the other 44 behave like LMXBs. 

Pietsch et al. (2005) conducted a survey of the X- 
ray population in all archival XMM-Newton observa- 
tions of M31. They detected 856 point sources, includ- 
ing 7 X-ray binaries and 9 candidate X-ray binaries. 
They also reported 27 X-ray source associations with 
globular clusters, with 10 further possible associations; 
these are also likely LMXBs. Of these 53 sources, 21 
are in our sample and 10 of those exhibited Type A 
PDS. Hence, 36 of our X-ray binaries are newly iden- 
tified in this work. 

In addition, luminosity variations greater than 5cr 
between observations were exhibited by 57 of the 
sources. Of these, 7 are transient, in that they vary in 
intensity by a factor of <;10, and are absent in at least 
one observation. The fraction of variable X-ray sources 
is 90% in our survey, higher than the 50% quoted by 
K02 for the field population; this is simply because the 
XMM-Newton observations were deeper than the ones 
performed by Chandra. Trudolyubov & Priedhorsky 
(2004) analysed XMM-Newton observations of M31 
X-ray sources in globular clusters, and found 80% to 
be variable. They suggested that the variability for 
LMXBs in globular clusters was higher than those 
in the field, comparing XMM-Newton observations of 
globular cluster X-ray sources with Chandra obser- 
vations of the general population. Our XMM-Newton 
observations suggest that the globular cluster X-ray 
sources are no more variable than the general popula- 
tion. 

Among our —50 likely LMXBs, we have observed 
1 Z-source, 2 dippers, 7 transients and 13 candidate 
black holes. Out of the -150 Galactic LMXBs, 6 or 
7 are Z-sources, —10 are dippers, —70 are transients 
and 37 are black hole candidates (e.g. Liu et al. 2001). 
Hence, our sample contains —half the number of Z- 
sources and dippers, and — one third of the transients 
expected from the Galactic population. However, we 
found a comparable number of black hole candidates. 
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Appendix A: Survey results in detail 

For each observation of every source, we recorded the 
fractional r.m.s variability, F val (Vaughan et al. 2003), 
PDS characteristics and SED characteristics. 

We obtained best fit parameters to the 64- and 
128- bin PDS using power law and broken power law 
models, and used F-testing to ascertain the probability 
that the broken power law is the better fit. The null 
hypothesis probability, Hq, is the probability that the 
deviation of observed data from the model is solely 
due to random statistical fluctuations. We rejected any 
fits where Hq <0.001, and considered any model with 
Ho >0.05 an acceptable fit. Power law fits to our 64 
and 128 bin PDS have 9 and 12 degrees of freedom, 
respectively, while broken power law fits have 8 and 
11 degrees of freedom for our 64 and 128 bin PDS. We 
present the % 2 corresponding to H = 0.001 and 0.05 
for 8, 9, 11 and 12 degrees of freedom in Table A.l. For 
example, we reject fits with x 2 /dof >26/8 and accept 
fits where x 2 /dof <16/8. 

We fitted each SED with the following sim- 
ple spectral models: blackbody (BB), bremsstrahlung 
(BR) and power law (PO), all suffering line-of-sight 
absorption. In some cases, a second component was re- 
quired; in this case, we used a blackbody + power law 
(2C) model. In addition, the supersoft sources were 
fitted with neutron star atmosphere (NSA) models; 
this was to represent a thermal model with a sharper 
drop off at low energies than a blackbody. 

We present the properties of each source in 
Table A. 2; we present the results of each observation 
on a different line. Numbers in parentheses indicate 
uncertainties in the last digit. Column 1 is the source 
name. Column 2 is the fractional r.m.s variability of 
the background subtracted lightcurve. Columns 3-7 
describe the fits to either the 64- or 128 bin PDS; 
we present results for the PDS with the lowest F-test 
probability, to give the most conservative detections 
of Type A variability. Columns 3 and 4 are the spec- 
tral index, 7 and x 2 /dof for the power law model. 
Similarly, columns 5, 6 and 7 are the break frequency, 
v c , spectral index for v > v c , (i, and x 2 /dof for the 

Table A.l. Limits in x 2 corresponding to probabilities of 
0.001 and 0.05, for a given number of degrees of freedom 



D.O.F x 2 {Ho = 0.001) x 2 (Ho = 0.05) 



broken power law model. Uncertainties in columns 3 
7 are at the 1 a level. The best fit model is indicated 
in boldface. Column 8 gives the F-test function, F, 
which is the probability that the improvement in fit- 
ting a broken power law to the PDS over a power law 
is due to chance. Column 9 is the PDS type. Columns 
10-14 describe the SED. Column 10 gives the best fit 
spectral model. Column 11 is the equivalent line-of- 
sight absorption, Ah, in units of 10 22 H atom cm~ 2 . 
Column 12 describes the shape of the SED, in terms 
of the temperature, kT, in kcV units for thermal mod- 
els (BB, BR or NSA), or the photon index, T, for the 
PO or 2C models. Column 13 gives the x 2 /dof for the 
best fit to the SED, while S14 gives the 0.3-10 keV 
luminosity, in units of 10 37 erg s _1 . Uncertainties in 
columns 10-12 and 14 are given at the 90% confidence 
level, and frozen parameters are indicated by "f" . 

The PDS are classified as follows. If F <0.001, the 
PDS is of Type A. If F > 0.001, but H <0.001 for the 
best power law fit and H >0.05 for the broken power 
law fit, then we classify the PDS as Type A*, and 
treat it as Type A. If, however, the PDS is of Type A 
but dominated by variability in the background, then 
we denote it as Type A* in Table A. 2, and treat it as 
unclassified. If a power law gives x 2 /dof SSI.0, and 7 
is >3cr above zero, then the PDS is of Type B. If the 
PDS docs not rise significantly above a power of 2, 
then we classify it as a flat PDS, indicated by — . Any 
PDS not classified as above is unclassified, indicated 
by a question mark. 



8 26 16 

9 28 17 

11 31 20 

12 33 21 
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Table A. 2. PDS and SED properties for each observation of every source. Observations 1-4 are described on separate 
lines. Numbers in parentheses indicate uncertainties in the last digit. Columns 1 and 2 give the source name and 
fractional r.m.s. variability. Columns 3-9 concern fits to the PDS; 3 and 4 give the spectral index and x 2 /dof for the 
best fit power law model to the PDS, while 5, 6 and 7 show the break frequency in mHz, spectral index for v > v c and 
X 2 /dof for the best broken power law fit to the PDS; column 8 is F-test probability that the improvement in fitting 
the broken power law over the power law is due to chance, and column nine gives the PDS Type. Uncertainties in fits 
to the PDS are given at the lcr level. Columns 10-14 describe fits to the SED; column 10 gives the best fit model (see 
text for abbreviations); 11 gives the equivalent Hydrogen absorption in the line of sight; 12 gives either the photon 
index, or temperature, appropriate to the model; 13 gives the xVdof for the best fit model and 14 gives the 0.3-10 
keV luminosity. Uncertainties in fitting the SED are at the 90% confidence level. ". . ." indicates that no good SED was 
obtained; either the source was not observed with the PN, or the SED had negative counts. 
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101/74 


24.6(17) 




0.07(14) 


0.13(2) 


12/9 


60(60) 


0.4(3) 


12/8 


NA 




_2 


PO 


0.49(5) 


1.48(9) 


158/149 


23.5(14) 


10 


0.13(5) 


0.12(2) 


18/12 


54(6) 


0.6(2) 


12.11 


0.04 




? 


PO 


0.24(4) 


2.71(11) 


70/49 


6.93(14) 




0.1(2) 


0.08(3) 


19/12 


36.8(2) 


0.20(10) 


20/11 


NA 




? 


PO 


0.23(5) 


2.8(2) 


16/20 


8.1(8) 




0.13(13) 


0.02(2) 


15/9 


57(2) 


0.27(11) 


10/8 


0.08 






PO 


0.26(5) 


3.1(3) 


45/65 


4.4(4) 




0.09(6) 


0.09(2) 


45/9 


35(6) 


0.28(6) 


13/8 


2 E- 


-3 


At 


PO 


0.29(2) 


3.14(11) 


124/89 


7.8(5) 



a iVH is in units of 10 H atom cm - ' 1 

^Either spectral index or temperature, depending on the mode; kT is in units of keV 
C L is in units of 10 37 erg s _1 

d These luminosities assume S6 is in M31, however it is a foreground star, hence its luminosity is reduced by ~6 orders 
of magnitude. 
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Table A. 2. continued 



Power Law Fit Broken Power Law Fit SED FIT 



S Fvelt 


7 


xl 






xl 


F 




Type Spec 7V H 


T/kT 


xl 


L 


11 (J.OOyO) 


n 41 791 


^0 /Q 




n ^9f4i 

\J . 'J Li I 


1Q/G 


0.04 




? 


PO 


0.09(3) 


1 4fl7Ql 


^1 /^Ll 

Ol/ Ol 


Q 4771 


0.03(24) 


0.10(3) 


20/9 


21(12) 


0.23(11) 


6/8 


0.07 




? 


PO 


0.10(5) 


1.33(13) 


28/21 


7.3(9) 


0.14(6) 


0.08(2) 


10/9 


34(12) 


0.21(10) 


6/8 


0.08 






PO 


0.12(2) 


1.56(8) 


86/74 


4.4(3) 


0.10(3) 


0.13(2) 


27/12 


22(4) 


0.22(4) 


11/11 


2 E- 


-3 


? 


PO 


14791 


1.59(6) 


167/135 


7.0(3) 


1 9 fl 047Q1 


n i n79l 


22/12 




n 4nCi i "i 


i n /i i 


3xl0~ 3 


? 


PO 


0.38(6) 


9 91 71 'il 


4^/43 

*LO j 4:0 


4 87^1 


0.17(11) 


0.09(4) 


9/9 


38(14) 


0.31(16) 


5/8 


0.02 




? 


PO 


0.38(11) 


2.2(3) 


11/13 


5.0(6) 


0.11(9) 


0.12(2) 


17/9 


30(70) 


0.22(5) 


14/8 


0.23 




? 


PO 


0.33(4) 


2.14(13) 


104/91 


4.6(3) 


0.15(4) 


0.17(2) 


32/9 


29(3) 


0.37(5) 


5/8 


9 T? 


A 




pn 




2.18(9) 


76/98 


6.0(3) 


i q n i 87i 91 

lO U. LoyLZ. J 


n nQ79l 


31/19 




n 4^7*41 


1 8 /I 1 


0.02 




? 


PO 


0.13(6) 


1 8791 


9^ /9/L 

ZO / Z1 


1 7791 


0.28(18) 


0.11(4) 


6/12 


50(20) 


0.3(3) 


11/11 


NA 




? 


PO 


0.19(15) 


2.6(8) 


3/4 


1.6(2) 


0.15(8) 


0.11(2) 


13/9 


36(8) 


0.34(10) 


2/8 


2 E- 


•4 


A 


PO 


0.13(3) 


1.85(14) 


62/71 


2.6(2) 


0.14(5) 


0.16(2) 


47/9 


29(3) 


0.37(5) 


8/8 


9 T? 


A 

-4 




JT W 


n i R731 


1.99(9) 


95/75 


3.3(4) 

\ / 


14 U.UU^ / I 


n 0^791 


R /Q 






7 /S 


0.40 




— 


PO 


0.19(2) 


1 73781 


8^ /81 
oo /or 


7 7741 


0.07(10) 


0.07(3) 


15/12 


50(20) 


0.3(3) 


13/11 


0.21 




? 


PO 


0.13(7) 


1.7(2) 


9/11 


5.7(7) 


0.11(2) 


0.13(2) 


23/12 


32(8) 


0.28(9) 


26/11 


NA 




? 


2C 


0.22(10) 


2.39(7) 


376/347 


24.6(12) 


0.099(14) 


0.18(2) 


15/12 


19(3) 


0.26(3) 


14/11 


n /in 




? 


9^ 


fl 93791 


2.33(6) 


373/375 


23.8(5) 


K n 1 o/i o\ 

LO U.IO^IOJ 


n 1 9791 

O. J.Z ^z ^ 


Q /Q 


9Q781 
zy yo J 


n 91 771 

U.Z± 1 t 1 


i n /s 


NA 




B 


PO 


0.32(9) 


9 7)731 


98 /on 


1 Q731 


0.1(4) 


0.11(4) 


7/9 


50(20) 


0.3(3) 


9/8 


NA 




? 


PO 


0.17(11) 


2.4(5) 


4/8 


1.2(3) 


0.17(16) 


0.08(2) 


18/12 


37(2) 


0.26(6) 


11/11 


0.02 




? 


PO 


0.36(8) 


2.3(2) 


43/43 


1.9(3) 


0.33(9) 


0.17(2) 


22/9 


29(4) 


0.33(5) 


9/8 


7 T? 


Q 
O 


? 


PO 


39781 


2.3(2) 


17/25 


1.4(2) 


1 fi R4781 


n 1 4791 

u. L^Lyz. f 


1 'i/Q 

LO / U 


32(6) 


0.32(7) 


8 /8 


0.02 




? 


PO 


0.2(2) 


9 1 771 


38/33 

OO / OO 


1 c i74l 


0.4(3) 


0.13(4) 


3/9 


28(14) 


0.23(12) 


5/8 


NA 




B 


PO 


0.2(2) 


1.8(7) 


4/8 


1.4(4) 


0.59(14) 


0.15(2) 


26/12 


51(6) 


0.63(17) 


24/11 


0.38 




? 


PO 


0.4(2) 


2.6(6) 


37/42 


1.9(5) 


0.30(11) 


0.15(2) 


33/9 


25(4) 


0.2(3) 


5/8 


1 Th 1 


A 


A 
/\ 


pn 


(1 9S7Q1 


1.8(2) 


35/33 


1.8(2) 

V / 


17 fl 1 47Q1 




14/Q 
11/ y 




n i R741 


1 1 /s 


0.23 




? 


PO 


0.10(4) 


1 ^(lA^ 
1 .OO y 14 J 


3^ /31 
oo /or 


3 4731 


0.02(76) 


0.03(3) 


13/12 


50(30) 


0.2(3) 


13/11 


0.78 






PO 


0.20(8) 


1.72(15) 


8/8 


3.9(6) 


0.13(9) 


0.10(2) 


15/9 


34(10) 


0.26(9) 


7/8 


0.02 




? 


PO 


0.14(3) 


1.69(12) 


62/64 


3.3(3) 


0.17(5) 


0.16(2) 


31/9 


29(4) 


0.33(5) 


6/8 


4 E- 


•4 


A 
/\ 


pn 


n i ^731 


1.67(9) 


58/61 


3.4(2) 


18 n on /'i o~\ 




1 1 / 1Z 


o^/'s^ 


n i of ^ 
u.iyi o i 


1 C /I 1 


0.24 




? 


PO 


0.122 


1 7/1771 




o.y 1^0 ) 


0.08(6) 


0.09(3) 


12/9 


26(10) 


0.20(7) 


10/8 


0.20 




? 


PO 


0.11(3) 


1.69(12) 


29/29 


8.8(8) 


0.09(4) 


0.12(2) 


24/12 


30(7) 


0.31(8) 


15/11 


0.03 




? 


PO 


0/114(9) 


1.74(6) 


189/160 


8.2(3) 


0.10(2) 


0.11(2) 


17/12 


32(6) 


0.23(5) 


19/11 


NA 




V 


PO 


0.12(2) 


1.79(5) 


146/173 


8.7(3) 


19 0.20(12) 


0.07(2) 


13/12 


32(14) 


0.14(8) 


15/11 


NA 




B 


PO 


0.157 


2.9(4) 


23/37 


3.3(4) 


0.1(4) 


0.11(4) 


15/9 


43.6(12) 


0.34(13) 


16/8 


NA 




? 


PO 


0.1 f 


2.3(2) 


10/9 


2.9(6) 


0.20(9) 


0.13(2) 


9/9 


22(8) 


0.22(5) 


6/8 


0.07 




B 


PO 


0.20(5) 


3.0(3) 


130/105 


5.1(5) 


0.11(4) 


0.16(2) 


40/9 


22(3) 


0.31(4) 


9/8 


6 E- 


-4 


A 


PO 


0.23(2) 


2.47(9) 


182/115 


9.7(5) 


20 0.07(7) 


0.01(3) 


6/9 


44(2) 


0.42(7) 


6/8 


0.63 






PO 


0.14(2) 


1.86(9) 


85/86 


7.8(5) 


0.10(7) 


0.14(3) 


8/9 


35(13) 


0.4(2) 


7/8 


0.23 




? 


PO 


0.13(4) 


1.83(14) 


47/33 


8.0(6) 


0.09(5) 


0.10(2) 


8/9 


35(10) 


0.26(9) 


6/8 


0.15 






PO 


0.14(2) 


1.84(7) 


142/127 


7.1(3) 


0.08(3) 


0.11(2) 


20/9 


38(5) 


0.30(6) 


9/8 


0.02 




? 


PO 


0.16(2) 


1.80(5) 


223/201 


9.1(3) 
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Table A. 2. continued 



Power Law Fit Broken Power Law Fit SED FIT 



S Fv&r 


7 


xl 




P 


xl 


F 




Type Spec Nu 


r/kT 


xl 


L 


21 27(17^ 

— ■ _L ' ' . _ I 111 I 


0.11(2) 


21/9 


1 QC31 


22(31 

V } % AM AM 1 « J 1 


7/8 


6 E- 


-3 




PO 


0.1 f 




29/18 


1 0C21 


0.21(17) 


0.06(4) 


5/9 


20(20) 


0.09(8) 


5/8 


NA 






PO 


0.6(2) 


2.6(4) 


9/5 


3.2(6) 


0.14(8) 


0.10(2) 


36/12 


48(7) 


0.58(17) 


10/11 


2 E- 


-4 


A t 

A 1 


PO 


0.28(4) 


2.0(13) 


62/71 


3.4(3) 


0.20(5) 


0.15(2) 


47/9 


37(4) 


0.44(6) 


11/8 


2 E- 


-3 


A* 


PO 


0.28(4) 


2.15(9) 


45/50 


2.2(3) 


99 n i ^f\} 


n n7(2l 


i n /q 




n i q(i ^1 

\j. l y ^ lo j 


1 9 /S 


NA 




B 


2C 


0.1 f 




77 /73 




0.16(14) 


0.16(4) 


11/9 


46(10) 


0.7(3) 


4/8 


7 E- 


-3 


? 


PO 


0.1 f 


1.4(2) 


9/15 


5.1(14) 


0.10(5) 


0.11(2) 


38/12 


26(6) 


0.27(6) 


24/11 


0.03 




*? 


2C 


0.1 f 


2.42(8) 


127/128 


6.6(12) 


0.16(3) 


0.13(2) 


53/9 


29(4) 


0.33(5) 


7/8 


3 E- 


-5 


At 


9P 


1 f 


2.40(6) 


139/114 


6.8(5) 




n 0^(91 


1 Q /I 9 


OU.OIO 1 




19/11 

J. £i j A. A. 


u.uo 




? 


PO 


0.06(3) 


9 1 A(ft 1 i 


43 l^tA 




10(10) 


0.22(3) 


38/8 


45(6) 


1.0(3) 


9/8 


9 E- 


-4 


A 












5(5) 


0.01(2) 


14/9 


33(8) 


0.24(5) 


4/8 


2 E- 


-3 


? 












0.9(16) 


0.04(2) 


20/12 


17(8) 


0.09(3) 


16/11 


0.12 




? 












94 n run nisi 


n n2(2l 


7 /Q 


fin(Q0l 




7 /8 


NA 






PO 


0.20(11) 


1 7(^ 


11/17 




0.31(19) 


0.06(3) 


13/12 


30(20) 


0.18(11) 


11/11 


0.15 




— 


PO 


0.15(14) 


1.5(4) 


18/15 


1.3(3) 


0.31(15) 


0.04(2) 


9/12 


36(17) 


0.12(8) 


9/11 


0.30 






PO 


0.32(8) 


1.9(3) 


46/38 


1.5(2) 


0.27(9) 


0.18(2) 


57/9 


31(3) 


0.45(5) 


10/8 


2 E- 


-4 




PO 




1.69(15) 


23/20 


1.11(2) 


9 c i f)9C4l 


n 04(91 


m/Q 


52(12) 


n 2^(1 61 

Ui^JI 1U 1 


fi /s 


n ofi 






2C 


0.152(5) 


9 1 1 


493/41IS 




0.02(3) 


0.06(3) 


10/12 


40(60) 


0.08(17) 


12/11 


NA 






2C 


0.134(11) 


2.02(6) 


189/204 


63(4) 


0.038(12) 


0.19(2) 


18/9 


0.02(4) 


0.32(5) 


7/8 


9 E- 


-3 


? 


2C 


0.162(12) 


2.23(3) 


715/666 


57(3) 


0.093(9) 


0.19(2) 


27/9 


34(4) 


0.94(7) 


21/8 


0.15 




B 


2C 


0.13(9) 


2.13(8) 


598/557 


51(2) 


2fi 1 (41 


-0 02(21 


11/12 


57.8(11) 


1 2(1 31 

\J • _|_ AM 1 _l_ « J 1 


10 /ll 


0.58 






PO 


0.19(9) 


2 7('3 1 l 


36/33 


1 5f3l 


0.3(8) 


0.15(4) 


9/9 


44(2) 


0.40(12) 


9/8 


0.52 




Y 


PO 


0.15 f 


3.1(7) 


25/23 


0.9(5) 


0.3(13) 


0.05(2) 


25/12 


25(9) 


0.18(6) 


14/11 


0.08 




Y 


PO 


0.07(7) 


2.5(7) 


41/51 


0.43(10) 


0.26(13) 


0.08(2) 


23/12 


35(7) 


0.24(6) 


6/11 


2 E- 


-4 


A 1 " 


PO 




2.9(3) 


62/64 


1.7(2) 


97 fl 9(^1 


0.15(2) 


14/9 


36(6) 




8 /8 


0.04 




? 


PO 


0.3 f 




42/30 


n 7(a) 

U. 1 (1) 


2(10) 


0.04(3) 


15/9 


57.1(8) 


0.3(2) 


12/8 


0.19 




*? 


PO 


0.1 f 


1.3(13) 


24/14 


0.4(3) 


1.3(7) 


0.20(2) 


40/12 


30(4) 


0.45(6) 


19/11 


5 E- 


-3 


A* 


PO 


0.6(4) 


1.5(5) 


40/40 


0.7(2) 


0.42(11) 


0.15(2) 


36/9 


20(4) 


0.26(4) 


7/8 


5 E- 


-4 


A 


PO 


0.40(11) 


1.80(15) 


31/38 


1.4(2) 


98 1 1 (41 

J L^4;^ 


n i s(9l 


^A /Q 

04:/ y 




n <^a( i nl 


fi /£. 


A P 


-4 


A T 


PO 


0.2 f 


1.6(0) 


7 17 




0.4(6) 


0.07(4) 


22/9 


57.1(2) 


0.1(2) 


24/8 


NA 




? 


PO 


0.2 f 


1.7(7) 


1/3 


0.4(2) 


0.9(3) 


0.21(2) 


59/12 


27(3) 


0.55(7) 


5/8 


3 E- 


-4 


A 


PO 


0.20(16) 


2.0(4) 


27/24 


0.4(3) 


0.79(17) 


0.21(2) 


77/9 


34(3) 


0.53(6) 


18/8 


9 E~ 


-4 


A 


PO 


0.2(2) 


1.5(3) 


15/8 


0.52(10) 


29 0.2(2) 


0.001(2) 


10/12 


46.5(9) 


0.06(8) 


10/11 


0.42 






NSA 


0.1 f 


0.068(10) 


32/25 


0.6(5) 


0.35(17) 


0.15(3) 


24/12 


29(3) 


0.43(7) 


10/11 


2 E- 


-3 


? 


NSA 


0.1 f 


0.09(4) 


6/6 


0.8(3) 


0.32(16) 


0.11(2) 


34/12 


46(5) 


0.50(11) 


22/11 


0.03 




? 


NSA 


O.lf 


0.065(7) 


55/48 


0.6(4) 


0.61(16) 


0.21(2) 


73/9 


38(3) 


0.60(7) 


18/8 


1 E- 


-3 


A* 


NSA 


0.1 f 


0.068(12) 


39/28 


0.3(2) 


30 0.09(5) 


0.12(2) 


15/9 


26(13) 


0.23(4) 


10/8 


0.09 




? 


PO 


0.1 f 


2.13(9) 


52/41 


2.7(3) 


0.10(14) 


0.06(3) 


10/12 


17(13) 


0.12(5) 


8/11 


0.13 






PO 


0.1 f 


2.22(14) 


18/12 


2.8(4) 


0.08(11) 


0.01(2) 


8/9 


70(14) 


0.3(4) 


8/8 


NA 






PO 


0.13(3) 


2.43(16) 


79/81 


3.0(2) 


0.12(4) 


0.14(2) 


21/9 


14(3) 


0.20(2) 


9/8 


9 E- 


-3 


? 


2C 


0.06(3) 


1.88(10) 


133/94 


3.3(2) 
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Table A. 2. continued 





Power Law Fit 


Broken Power Law Fit 










SED FIT 






S Fvar 


7 


xl 


Vc 


P 


xl 


F 


Type Spec 


iV H 


T/kT 


xl 


L 


31 0.09(5) 


U.14U(z) 


OA /O 


48(7) 


n p^i (a 
U.o±(±o ) 


la/o 


0.20 


? 


2C 


0.25(8) 


2.1(2) 


86/90 


6.9(7) 


0.0b(o) 


0.10(3) 


26/9 


43(8) 


0.52(17) 


12/8 


n no 
U.U2 


? 


PO 


0.1o(3 J 


2.U(1U ) 


OD/OO 


IK 7/1 1 \ 

lo./(ll) 


0.06(2) 


0.11(2) 


26/12 


27(6) 


0.28(6) 


11/11 


2 E— 3 


? 


2C 


0.206(9) 


2.54(6) 


446/391 


23.9(10) 


0.07(2) 


0.15(2) 


39/9 


28(4) 


0.34(5) 


6/8 


2 E— 4 






0.20(2) 


2.40(5) 


353/314 


16.5(4) 


32 0.34(12) 


U. 10(Z ) 


OA /Q 




U.<±<±( 1U ) 


14/ o 


0.05 


? 


PO 


0.43(0.16) 


1.7(3) 


40/34 


3.6(5) 


0.18(10) 


0.13(3) 


14/12 


31(13) 


0.29(15) 


13/11 


0.DO 


? 


PO 


0.OD(lD J 


1.6(2) 


11 / 1 A 
11/ 10 


0.0(9) 


0.11(6) 


0.11(2) 


8/9 


17(2) 


0.18(3) 


4/8 


0.02 


B 


PO 


4(6) 


1.84(13) 


148/149 


9 8(7) 


0.04(12) 


0.07(2) 


11/12 


19(6) 


1.00(3) 


36/11 


NA 




PO 


0.42(3) 


1.95(8) 


120/106 


9 4(5) 


33 0.02(13) 


U.UZ \Z ) 


i a / ~i o 
lo/ 1^ 


Of? Q(E.\ 
OD.O^O ) 


U.Uo(o ) 


1 Q /I 1 

10/ 11 


NA 


~ 


PO 


0.09(2) 


1.87(7) 


146/153 


11.8(7) 


0.07(9) 


0.1(3) 


12/9 


14(9) 


0.23(5) 


11/8 


0.53 


? 


PO 


0.27(5) 


2.04(17) 


43/41 


8.4(8) 


0.08(4) 


0.16(2) 


8/9 


18(6) 


0.23(5) 


9/8 


NA 


B 


PO 


0.09(2) 


1.95(6) 


213/182 


10.7(5) 


0.10(2) 


0.23(2) 


48/9 


31(3) 


0.49(5) 


12/8 


1 E— 3 


A 


PO 


0.14(2) 


1 97f4) 


350/309 


10.7(3) 


34 0.05(6) 


0.02(2) 


5/9 


50(30) 


0.10(16) 


4/8 


0.47 




PO 


0.20(2) 


2.03(7) 


112/124 


15.3(6) 


0.04(11) 


0.13(4) 


18/9 


44(10) 


0.61(13) 


7/8 


6 E-3 


? 


PO 


0.18(4) 


1.96(12) 


27/35 


10.4(8) 


0.06(4) 


0.14(2) 


6/9 


27(7) 


0.26(7) 


7/8 


NA 


B 


PO 


0.16(2) 


1.78(5) 


172/165 


13.6(5) 


0.11(3) 


0.09(2) 


19/9 


39(6) 


0.27(7) 


5/8 


2 E-3 


? 


PO 


0.22(2) 


2.20(5) 


132/132 


7.0(3) 


35 0.08(5) 


0.16(2) 


14/9 


12(4) 


0.23(4) 


11/8 


0.24 


? 


NSA 


0.18(2) 


0.017(1) 


102/80 


38(3) 


0.12(8) 


0.12(3) 


23/9 


33(11) 


0.40(18) 


10/8 


0.02 


? 


NSA 


0.17(4) 


0.017(1) 


34/28 


26(6) 


0.10(6) 


0.07(2) 


17/9 


3(50) 


0.02(2) 


17/8 


NA 




NSA 


0.16(2) 


0.0164(4) 


86/76 


32(2) 


0.08(4) 


0.07(2) 


25/12 


27(2) 


0.14(3) 


18/11 


0.06 


? 


NSA e 


0.16(2) 


0.018(2) 


60/77 


24(7) 



36 ... 



0.3(6) 


0.11(2) 


13/9 


33(2) 


0.25(3) 


13/8 


NA 




I s 


PO 


0.27(7) 


1.5(2) 


46/29 


4.1(5) 


37 0.06(6) 


0.06(2) 


25/9 


50(12) 


0.29(15) 


20/8 


0.23 




? 


PO 


0.15(2) 


1.88(8) 


128/119 


9.3(6) 


0.02(19) 


0.02(3) 


16/12 


38(2) 


0.10(10) 


16/11 


0.43 






PO 


0.11(3) 


1.84(12) 


48/35 


8.9(7) 


0.12(4) 


0.131(14) 


34/12 


34(4) 


0.34(6) 


18/11 


0.01 




A* 


PO 


0.17(2) 


1.89(6) 


233/195 


8.6(3) 


0.08(3) 


0.21(2) 


51/9 


35(3) 


0.50(6) 


7/8 


2 E- 


-4 


A 


PO 


0.15(2) 


1.84(5) 


204/186 


8.7(3) 


38 0.19(6) 


0.22(2) 


19/9 


27(4) 


0.44(7) 


9/8 


0.02 




? 


PO 


0.15(4) 


1.72(17) 


109/85 


4.8(5) 


0.2(2) 


0.08(3) 


31/12 


44(8) 


0.54(19) 


16/11 


8 E- 


-3 


? 


PO 


0.13(13) 


2.2(5) 


4/5 


2.1(5) 


0.22(18) 


0.10(2) 


30/12 


44(5) 


0.44(11) 


14/11 


0.06 




? 


PO 


0.26(9) 


3.2(7) 


67/83 


2.7(4) 


0.24(5) 


0.18(2) 


50/9 


35(3) 


0.47(6) 


7/8 


2 E- 


-4 


A 


PO 


0.17(4) 


2.08(15) 


143/87 


2.8(3) 


39 0.2(3) 


0.14(2) 


29/12 


39(6) 


0.43(11) 


20/11 


0.04 




? 


PO 


0.1 f 


1.7(4) 


12/25 


0.6(3) 


0.09(51) 


0.13(3) 


18/9 


34(10) 


0.41(17) 


8/8 


0.02 




? 


PO 


0.1 f 


1.9(3) 


7/7 


1.2(5) 


0.1(6) 


0.07(2) 


12/9 


32(2) 


0.15(5) 


13/8 


NA 






PO 


0.10(6) 


1.80(16) 


24/49 


1.2(2) 


0.1(3) 


0.12(2) 


38/12 


31(5) 


0.32(6) 


6/11 


9E- 


5 


A* 


PO 


0.08(4) 


1.7(2) 


44/46 


0.97(14) 


40 0.15(6) 


0.12(2) 


18/9 


37(7) 


0.37(10) 


8/8 


0.01 




? 


PO 


0.15(3) 


2.02(13) 


60/45 


4.2(3) 


0.16(8) 


0.06(3) 


9/9 


51(16) 


0.3(2) 


7/8 


0.16 




? 


PO 


0.14(7) 


2.1(3) 


9/12 


4.3(5) 


0.09(6) 


0.18(2) 


56/9 


41(3) 


0.60(6) 


11/8 


4 E- 


-4 


A 


PO 


0.14(2) 


1.88(9) 


93/99 


4.9(3) 


0.11(5) 


0.06(2) 


14/12 


6.8(3) 


0.05(2) 


15/11 


NA 






PO 


0.15(2) 


2.09(8) 


85/92 


4.3(2) 



e Additional emission feature required, well fitted by 0.541(6) keV Gaussian with a = 0.036(8) 
/ Only PN data were available, not MOS 
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Table A. 2. continued. 



Power Law Fit Broken Power Law Fit SED FIT 



s 


Fv&T 


7 


xl 




P 


xl 


F 




Type Spec Nn 


T/kT 


xl 


L 


41 










O 33f 1 1 ^ 

U>OOI 11 I 


11/11 


n 1 s 




? 


PO 


0.06(6) 


L.I (O) 


97/90 






0.2(5) 


0.10(4) 


5/9 


37(16) 


0.26(15) 


3/8 


0.06 




? 


PO 


0.1 f 


2.1(5) 


5/4 


0.8(3) 




0.38(18) 


0.13(2) 


22/9 


42(5) 


0.39(9) 


17/8 


0.17 




? 


PO 


0.1 f 


2.3(4) 


23/22 


0.7(3) 




0.38(11) 


0.18(2) 


43/9 


26(3) 


0.36(4) 


6/8 


9 E- 


-5 


A 

A 


Tie \ 

PO 


0.09(5) 


1.71(15) 


51/43 


0.97(14) 


42 


0.13(6) 


0.09(2) 


25/9 


30(7) 


0.28(7) 


9/8 


5 E- 


-3 


? 


PO 


0.18(3) 


2.18(14) 


53/52 


4.4(3) 




0.10(11) 


0.05(4) 


18/9 


44(2) 


0.29(12) 


12/8 


0.10 




? 


PO 


0.20 


2.02(13) 


41/28 


5.0(7) 




0.04(17) 


0.08(2) 


8/9 


60(40) 


0.2(3) 


7/8 


0.3 






PO 


0.12(2) 


1.87(9) 


121/111 


4.9(3) 




0.12(3) 


0.17(2) 


34/9 


30(4) 


0.36(5) 


7/8 


6 E- 


-4 


At 


PO 


0.16(2) 


2.02(7) 


131/110 


5.2(2) 


43 0.09(10) 


0.06(2) 


7/9 


18(13) 


0.1(4) 


5/8 


0.13 






PO 


0.13(3) 


1.56(10) 


58/60 


4.8(4) 




0.02(40) 


0.18(3) 


19/9 


46(8) 


0.81(3) 


4/8 


4 E- 


-4 


A 


PO 


0.17(8) 


1.7(2) 


12/18 


5.1(6) 




0.10(6) 


0.11(2) 


39/12 


31(5) 


0.37(7) 


8/11 


4 E- 


-5 


A 


PO 


0.13(2) 


1.58(8) 


114/126 


5.5(3) 




0.14(3) 


0.20(2) 


47/9 


34(4) 


0.49(2) 


20/8 


0.01 




? 


PO 


0.14(2) 


1.58(6) 


166/138 


5.8(2) 



44 ... 



0.21(14) 


0.04(2) 


11/9 


14(8) 


0.05(2) 


12/8 


NA 




_/ 


PO 


0.12(5) 


0.98(10) 


47/48 


6.5(5) 


45 1.8(5) 


0.20(2) 


34/9 


45(5) 


0.68(17) 


12/8 


5 E- 


-3 


A* 


PO 


0.9(9) 


0.9(10) 


18/12 


0.5(3) 


0.5(12) 


0.08(3) 


15/9 


37(2) 


0.24(9) 


14/8 


0.45 




? 


PO 


0.1 f 


1.3(9) 


10/6 


0.4(3) 


1.0(4) 


0.08(2) 


4/9 


9(7) 


0.10(4) 


5/8 


NA 




B 


PO 


0.6(6) 


1.5(9) 


16/22 


0.5(2) 


3.5(6) 


0.23(2) 


91/9 


30(2) 


0.53(4) 


5/8 


3 E- 


-6 


A 


PO 


3(2) 


3(2) 


2/6 


0.4(2) 


46 0.3(2) 


0.09(2) 


4/12 


50(5) 


0.10(2) 


4/11 


NA 




B 


PO 


0.1 f 


1.9(4) 


9/7 


1.9(7) 


0.5(5) 


0.15(4) 


14/12 


46(14) 


0.3(2) 


17/11 


NA 




? 


PO 


0.1 


2.2(5) 


4/5 


0.6(3) 


0.2(2) 


0.03(2) 


14/12 


29(21) 


0.08(6) 


13/11 


0.41 




? 


2C 


0.1 f 


2.5(2) 


35/24 


2.32(10) 


0.55(10) 


0.22(2) 


41/9 


25(3) 


0.39(4) 


10/8 


1 E- 


-3 


A* 


PO 


0.1 f 


1.83(17) 


17/13 


0.59(10) 


47 0.13(6) 


0.14(2) 


15/9 


37(7) 


0.34(10) 


11/8 


0.14 




? 


PO 


0.10(3) 


2.28(15) 


65/54 


3.5(3) 


0.09(12) 


0.07(4) 


12/9 


19(20) 


0.05(10) 


13/8 


NA 






PO 


0.09(5) 


1.9(2) 


20/18 


4.6(6) 


0.13(6) 


0.15(2) 


42/9 


41(4) 


0.59(9) 


5/8 


6 E- 


-5 


A 


PO 


0.12(2) 


2.01(11) 


110/97 


3.4(2) 


0.11(2) 


0.21(2) 


62/ 9 


21(3) 


0.38(4) 


11/8 


3 E- 


-4 


A 


PO 


0.15(2) 


1.89(4) 


265/231 


11.9(3) 


48 0.19(9) 


0.09(2) 


22/9 


32(7) 


0.28(7) 


8/8 


7 E- 


-3 


? 


PO 


0.09(4) 


2.2(2) 


29/24 


1.9(3) 


0.5(2) 


0.13(3) 


23/9 


48(8) 


0.7(3) 


10/8 


0.01 




? 


PO 


0.1 f 


2.1(3) 


12/10 


1.1(3) 


0.18(10) 


0.12(2) 


63/9 


42(2) 


0.51(5) 


21/8 


4 E- 


-3 


? 


PO 


0.11(4) 


2.3(3) 


53/50 


1.6(2) 


0.36(9) 


0.19(2) 


42/9 


22(3) 


0.34(4) 


6/8 


1 E- 


-4 


A 


PO 


0.14(4) 


2.2(2) 


41/35 


1.6(2) 


49 0.2(2) 


0.08(2) 


16/9 


35(10) 


0.24(9) 


8/8 


0.03 




? 












0.17(9) 


0.15(4) 


12/8 


48(10) 


0.6(3) 


5/8 


0.01 




? 


PO 


0.12(8) 


1.9(3) 


11/13 


3.9(6) 


0.20(11) 


0.15(2) 


17/9 


41(7) 


0.34(9) 


26/8 


NA 




? 












0.17(5) 


0.19(2) 


63/9 


37(3) 


0.53(6) 


12/8 


4 E- 


-4 


A 


PO 


0.13(4) 


1.87(13) 


52/37 


2.4(2) 


50 0.1(10) 


0.07(2) 


33/12 


27(7) 


0.27(6) 


13/11 


1 E- 


-3 


A* 


PO 


0.1 f 


2.2(4) 


8/14 


0.9(3) 


0.1(9) 


0.18(3) 


14/9 


20(8) 


0.31(9) 


11/8 


0.21 




? 


PO 


0.3(2) 


3.1(9) 


3/7 


1.2(3) 


0.33(10) 


0.10(2) 


5/9 


3(40) 


0.1(2) 


5/8 


NA 




B 


PO 


0.07(4) 


1.78(14) 


39/40 


1.07(10) 
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Table A. 2. continued. 



Power Law Fit Broken Power Law Fit SED FIT 



S Fvelt 


7 


xl 


v c 


P 


xl 


F 




Type Spec 
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r/kT 


xl 


L 
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08(91 
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11/11 


0.22 




? 


NSA 


0.1 f 
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U . UOl u / 


29 /25 




3(2) 


0.22(3) 


11/12 


24(6) 


0.47(12) 


4/11 


6 E- 


•4 


A t 

A T 


NSA 


0.1 f 


0.08(5) 


3/5 


0.8(8) 


1.4(12) 


0.14(2) 


23/9 


48(6) 


0.58(15) 


4/8 


2 E- 


•4 


A 

A 


PO 


0.1 f 


4(4) 


30/23 


0.4(4) 


0.0.09(3) 


0.12(2) 


19/9 


37(6) 


0.30 (7) 


10/8 


0.03 




? 


2C 


0.08(2) 


2.30(11) 


150/147 


5.5(3) 


59 ()AQ(~\9\ 


0.19(2) 


32/9 




0.47(8) 


13 /8 


0.01 




A* 


PO 


0.1 f 


2.2(6) 


18/12 


u. / yo) 


0.4(3) 


0.07(3) 


11/12 


0.56(12) 


0.5(3) 


10/11 


0.73 






PO 


0.1 f 


0.9(4) 


6/6 


1.2(6) 


0.41(13) 


0.10(2) 


8/12 


16(7) 


0.15(4) 


7/11 


0.29 




B 


PO 


0.1 f 


1.2(2) 


47/31 


1.1(3) 


0.45(10) 


0.15(2) 


27/9 


23(4) 


0.25(4) 


10/8 


6 E- 


•3 


? 


PO 


0.1 f 


1.08(15) 


29/20 


0.93(17) 


JU 1.41 O i 


()P,(9) 


9 /9 


46(15) 


D 2<Ul 4"! 


6 /8 


0.09 




? 


PO 


0.1 f 


3.5(8) 


21/19 


0.21(14) 


11(14) 


0.15(3) 


33/9 


39(7) 


0.61(17) 


8/8 


9 E- 


•4 


A t 

A T 


NSA 


0.1 f 


0.05(5) 


12/8 


0.2(2) 


2.0(14) 


0.10(2) 


10/12 


27(2) 


0.18(4) 


6/8 


0.01 




? 


PO 


0.1 f 


3.6(13 


19/19 


0.16(12) 


0.18(4) 


0.21(2) 


35/9 


24(3) 


0.37(4) 


11/8 


3 E- 


•3 


A * 


PO 
I w 


0.17(3) 


1.83(8) 


88/85 


4.2(2) 




n 1 c \(9\ 


fi /Q 


yui ou i 




o/ o 


u.uo 




? 


PO 


n nsnn a\ 




1 1 8 /I 9Q 




0.01(29) 


0.09(4) 


15/9 


28(14) 


0.24(12) 


12/8 


0.14 




? 


PO 


0.11(3) 


1.71(11) 


37/43 


13.2(9) 


0.07(4) 


0.03(2) 


13/9 


12(2) 


0.05(3) 


12/8 


0.47 






PO 


0.10(2) 


1.79(4) 


204/203 


11.4(4) 


0.08(12) 


0.19(2) 


46/9 


38(3) 


0.52(6) 


18/8 


7 E- 


•3 


? 


PO 


0.11(2) 


1.72(2) 


216/200 


10.3(3) 




0.11(2) 


9 /9 


1 A(9\ 




i n/s 


NA 




B 


PO 


n in(5l 


n (\a(K\ 

.u^yo j 




±4;. y 1 -L -L J 


0.12(7) 


0.16(3) 


20/12 


37(2) 


0.49(8) 


10/11 


6 E- 


•3 


? 


2C 


0.2(2) 


2.1(9) 


22/23 


21(6) 


0.07(6) 


0.10(2) 


6/9 


27(8) 


0.21(7) 


2/8 


0.01 




? 


2C 


0.1 f 


2.4(3) 


156/169 


22(3) 


0.45(4) 


0.17(2) 


34/12 


22(4) 


0.29(4) 


21/11 


0.02 




? 


PO 


0.05(4) 


0.52(7) 


89/72 


12.1(8) 


56 12dl1 


05C2) 


7/9 


33(16) 


0.14(5) 


6/8 


0.15 




? 


NSA 


0.24(5) 


0176('8~) 


51/49 


36(15) 


0.3(7) 


0.06(4) 


25/12 


20(20) 


0.15(11) 


22/11 


0.71 




? 


NSA 


0.3 f 


0.029(9) 


7/11 


3(2) 


0.9(5) 


0.08(2) 


20/12 


36(2) 


0.30(5) 


6/11 


5 E- 


•4 


A t 

A T 


NSA 


0.3 f 


0.027(11) 


18/14 


1.5(10) 


0.32(6) 


0.20(2) 


32/8 


25(3) 


0.35(4) 


7/8 


6 E- 


•4 
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? 
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^iUyo j 
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0.18(3) 
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23(7) 


0.44(12) 


16/11 


8 E- 


•3 


? 


NSA 


0.1 f 


0.030(4) 


4/10 


0.7(5) 


1.1(6) 


0.12(2) 


44(9) 


41(5) 


0.52(9) 


11/8 


1 E- 


•3 


A t 

A T 












0.45(11) 


0.17(2) 


68/9 


36(3) 


0.50(6) 


13/8 


4 E- 


•4 




PO 


0.05(5) 


2.3(3) 


36/41 


0.62(14) 


OO U.O-L I TO J 




i n /q 


35(13) 


n onfl n\ 

KJ ■ £t\J 1 1U 1 


10 /8 


0.22 




? 


BB 


n i f 

U. -L 1 




L ±a j OO 




0.3(3) 


0.18(4) 


10/9 


39(9) 


0.47(15) 


8/9 


0.86 




? 


BB 


0.1 f 


0.12(3) 


8/8 


0.8(6) 


0.2(3) 


0.07(2) 


11/9 


24(2) 


0.17(4) 


4/8 


3 E- 


•3 


? 


BB 


0.21(11) 


0.17(2) 


75/86 


2.1(4) 


0.32(8) 


0.20(2) 


88/12 


33(3) 


0.56(6) 


14/11 
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Power Law Fit Broken Power Law Fit SED FIT 

7 xl P Xl F Type Spec N H T/kT X l L 



Table A. 2. continued. 



61 0.23(9) 0.16(2) 9/9 27(6) 0.31(6) 4/8 
0.27(14) 0.12(3) 24/9 42(7) 0.58(17) 9/8 
0.23(10) 0.18(2) 31/9 31(4) 0.44(7) 5/8 
0.13(7) 0.12(2) 35/9 35(5) 0.32(6) 7/8 

62 0.09(7) 0.07(2) 15/12 38(11) 0.25(10) 6/11 
0.12(10) 0.16(4) 20/9 42(6) 0.60(17) 10/8 
0.08(7) 0.12(2) 14/9 57(5) 0.6(3) 21/8 
0.14(3) 0.18(2) 42/9 29(3) 0.40(4) 6/8 

63 ... ... 

0.26(16) 0.22(3) 23/12 31(7) 0.63(15) 6/11 
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0.21(2) 0.17(6) 30/12 25(4) 0.25(4) 8.3/11 2 E-4 A PO 0.23(3) 2.84(16) 66/65 6.4(5) 



